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INTRODUCTION 

1.1  Shock  Boundary-Layer  Interaction 
On  the  surface  of  high  speed  vehicles  in  transonic, supersonic  and  hypersonic 
flow,  shocks  generated  by  the  vehicle  geometry  can  interact  with  the  surface 
boundary  layers.  The  interaction  can  cause  a  number  of  problems  when 
designing  the  vehicles,  such  as  local  heating  (Fig. la)  and  boundary-layer 
separation.  The  local  heating  problem  sometimes  requires  the  vehicle  to 
employ  heat-resisting  materials  over  a  large  area  of  the  surface. 

These  problems  arise  because  high  flight  speeds  generate  high  temperatures 
in  the  boundary  layer  and  high  pressures  on  the  vehicle  surface.  Shock 
wave  boundary  layer  interaction  will  usually  increase  the  local  pressure 
still  further  and  this  often  implies  a  corresponding  increase  in  local  heat 
transfer  as  shown  in  Fig. lb,  taken  from  Neumann  (Ref. 19).  If  the  inter¬ 
action  is  strong  enough  to  cause  separation  then  the  whole  aerodynamic 
performance  may  be  affected  with  phenomena  such  as  buffet,  wing-rock  and 
wing-drop  limiting  manoeuvrability.  These  limitations  are  often  most 
severe  at  transonic  speeds  and  an  example  of  shock-induced  separation  is 
shown  in  Fig. 2  taken  from  Kubota  et  al.(Ref.2). 

In  the  case  of  two-dimensional  shock  boundary-layer  interaction  (Fig. 3, from 
Kuchemann  et  al..  Ref. 3),  theoretical  and  experimental  investigations  have 
had  some  success  in  confirming  the  general  flow  features  of  the  interaction. 
Fig. 4  (from  Horstman  et  al..  Refs.  4  and  5)  shows  some  examples  in  which 
there  is  good  agreement  between  the  experimental  results  and  theoretical 
predictions.  However,  for  three-dimensional  interaction  which  is  the  more 
practical  case,  there  are  still  a  number  of  unresolved  problems  because 


of  the  complicated  flow  structure  and  the  experimental  difficulties 
involved  in  the  full  resolution  of  the  physical  phenomena. 

Glancing  interaction  is  one  of  the  most  typical  three-dimensional  shock 
boundary- layer  interactions.  This  consists  of  an  oblique  shock  glancing 
across  a  boundary  layer  growing  along  a  wall  placed  at  right  angles  to  the 
shock-wave  surface.  In  most  of  the  experimental  investigations,  a 
turbulent  boundary  layer  growing  along  a  flat  plate  (called  the  side  wall 
in  this  paper)  is  employed  as  the  test  boundary  layer,  and  the  oblique 
shock  is  generated  by  a  wedge  mounted  normally  on  the  plate.  The  test 
arrangement  is  shown  in  Fig. 5  together  with  typical  patterns  of  flow 
properties  over  the  interaction  region. 

In  the  interaction  region  the  static  pressure  and  heat  transfer  distributions 
on  the  side  wall  are  very  similar  (Fig. 5).  For  a  weak  shock,  the  pressure 
and  heat  transfer  are  gradually  decreasing  with  increasing  distance  from  the 
shock  generator.  As  the  shock  strength  (pressure  ratio  across  the  shock)  is 
increased,  the  distributions  show  two  peaks,  one  near  the  shock  and  the 
other  near  the  shock  generator.  The  side-wall  boundary  layer  is  heavily 
skewed  near  the  shock,  because  the  surface  streamlines  are  deflected  much 
more  than  the  streamlines  at  the  edge  of  the  boundary  layer  which  only  turn 
to  lie  parallel  to  the  shock  generator  surface.  With  increasing  shock 
strength,  the  surface  flow  deflection  angle  increases  and  eventually 
exceeds  the  shock  angle.  As  a  result,  the  highly  deflected  surface-flow 
forms  a  particular  line  on  the  side  wall  upstream  of  the  shock  at  about  the 
shock-wave  angle.  In  front  of  this  particular  line,  the  surface  stream 
lines  tend  to  coalesce  with  it.  However,  behind  the  line,  the  surface 

stream  lines  run  nearly  parallel  to  it  or  they  approach  it  at  a 


small  angle.  This  coalescence  of  the  upstream  surface-stream  lines 
suggests  a  three-dimensional  flow  separation,  similar  to  Maskell's 
'ordinary  separation'  (Ref. 6).  As  shown  in  Fig. 6,  in  Maskell's 
definition  of  the  ordinary  separation,  the  separation  line  is  an  envelope 
of  surface  stream  lines  (or  surface  friction  lines)  converging  asymptotically 
from  each  side.  However,  the  surface  flow  found  in  this  interaction  is 
typical  of  what  most  investigators  regard  as  ordinary  separation  (refer 
to  the  oil-flow  picture  in  Fig. 5). 

One  of  the  earliest  experimental  studies  was  performed  by  Stanbrook  (Ref. 7). 
He  proposed  that  a  pressure  ratio  across  the  shock  of  1.5  was  sufficient  to 
induce  'separation'  at  Mach  numbers  from  1.6  to  2.0.  McCabe  (Ref. 8) 
developed  an  approximate  theory  for  the  surface-flow  deflection  in  a  quasi- 
two-dimensional  flow  field.  This  quasi-two-dimensional  flow  field  model 
is  an  idea  based  on  a  two-dimensional  interaction  caused  by  a  normal  shock 
(Fig. 7a).  It  is  reasonable  to  imagine  that  a  sectional  view  of  the  glancing 
interaction  region  normal  to  the  shock  is  similar  to  this  two-dimensional 
interaction.  It  may  also  be  proper  to  assume  that  all  flow  characteristics 
are  locally  regarded  as  constant  along  lines  parallel  to  the  shock  (Fig. 7b). 
Thus,  the  quasi-two-dimensional  flow  field  model  is  introduced  to  give  the 
idealized  surface-flow  patterns  illustrated  in  Fig. 7c.  McCabe  used  this 
flow  field  for  his  theory  and  further  assumed  that  all  vorticity  in  the 
side-wall  boundary  layer  is  convected  with  the  free  stream  velocity. 

McCabe's  theoretical  surface-flow  deflection  angle  generally  agrees  well 
with  the  experimental  maximum  surface-flow  angle,  when  the  shock  is  not  too 
strong.  As  an  extension  of  this  theory,  he  also  tried  to  predict  the  shock 
generator  angle  to  induce  'separation',  by  introducing  a  criterion  derived 


from  Che  geometrical  considerations  of  the  quasi-two-dimensional  flow  field. 

He  predicted  that  the  'separation*  is  incipient  when  the  deflected  surface- 
flow  becomes  aligned  with  the  shock  (Fig.7c-2).  This  prediction  gives 
slightly  higher  shock  generator  angles  than  his  experimental  results  at 
Mach  numbers  of  2  and  3.  Korkegi  (Ref. 9)  generalized  McCabe's  theoretical 
prediction  by  using  a  first  order  approximation  to  give  a  term  being  roughly 
proportional  to  the  shock  strength.  He  proposed  that  'incipient  separation' 
takes  place  when  the  product  of  the  free  stream  Mach  number  and  the  shock 
generator  angle  reaches  0.364  radians  (this  corresponds  to  a  pressure  ratio 
across  the  shock  of  1.63).  He  also  suggested  that,  to  obtain  better 
agreement  with  experimental  results,  the  product  of  the  Mach  number  and 
the  shock  generator  angle  should  be  approximately  0.3  (which  corresponds 
to  a  pressure  ratio  of  1.5).  These  mathematical  and  empirical  criteria 
are  shown  in  Fig. 9,  compared  with  some  experimental  incipient  separation 
data  defined  according  to  McCabe's  criterion.  As  shown  in  the  figure, 
these  criteria  roughly  coincide  with  the  experimental  results.  In  other 
words,  these  criteria  provide  a  good  prediction  of  the  shock  generator 
angle  at  which  the  surface-flow  deflection  becomes  equal  to  the  shock  angle. 
However,  they  do  not  provide  any  clarification  of  'incipient  separation' 
itself. 

Recent  experimental  investigations  have  supplied  a  great  deal  of  flow  field 
data  for  the  ineraction  region.  Many  interpretations  of  the  flow  field  data 
have  brought  diverse  explanations  of  the  interaction  phenomena.  For  instance, 
the  'corner 'flow  in  the  interaction  region  is  not  quasi-two  dimensional 
but  'conical'  as  shown  in  Fig. 8,  so  McCabe’s  criterion  may  not  be  applied 
to  define  'incipient  separation'.  This  is  because  the  'conical' 
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flow  field  allows  an  attached  convergent  surface-flow  even  when  the  surface 
flow  deflection  exceeds  the  shock  angle  (Fig. 9c).  The  occurrence  of 
'separation'  is  just  one  of  the  points  of  controversy. 

Oskam  et  al.  (Refs.  10-12)  discussed  in  some  detail  the  uncertainty  in 
McCabe's  criterion,  mostly  on  the  interpretation  of  their  own  oil-flow 
pictures  accompanied  by  viscous  layer  surveys.  They  concluded  that,  even 
when  the  surface-flow  deflection  exceeded  the  shock  angle,  the  interaction 
region  still  kept  an  attached  flow  pattern.  Although  the  oil  streak  lines 
were  converging  and  tended  to  coalesce  upstream  of  the  shock,  the  lines 
did  not  converge  into  a  single  line  and  the  region  in  which  the  oil 
coalescence  occurred  was  growing  in  size  in  a  direction  along  the  shock. 
This  was  consistent  with  the  idea  of  an  interaction  region  increasing  its 
size  with  distance  along  the  shock. 

Peake  (Refs. 13  and  14)  presented  a  different  view  based  on  his  detailed 
measurements  of  wall  surface  flow  patterns,  surface  static  pressures  and 
viscous  layer  properties.  In  particular  he  noticed  the  coalescence  of  oil- 
flow  in  a  region  where  the  surface  static  pressure  peaked.  He  therefore 
suggested  that  the  'three-dimensional  separation  line'  produced  a  free 
shear  layer  which  rolled  up  into  a  flattened  vortex  within  the  depth  of 
the  side-wall  boundary  layer  and  that  there  appeared  to  be  no  sudden 
eruption  of  vortical  fluid  from  the  surface.  He  also  remarked  that  the 
development  of  the  viscous  layers  with  increasing  shock  generator  angle 
was  so  gradual  that  it  was  difficult  to  pin  down  the  exact  shock  generator 
angle  needed  to  induce  'incipient  separation’  without  adopting  McCabe's 
criterion. 


Neumann,  Hayes  and  Token  (Refs. 15-20)  mainly  concentrated  on  investigations 
of  the  heat  transfer,  vhich  is  very  important  for  direct  practical 
applications.  According  to  their  studies,  the  heat  transfer  distributions 
in  a  direction  normal  to  the  free  stream  show  a  good  correlation  with  the 
surface  static  pressure  distributions,  as  shown  in  Fig. 5b.  These  distributions 
also  have  a  very  high  'peak'  near  the  shock  generator  root,  and  sometimes 
show  a  small  'peak'  close  to  the  surface-flow  coalescence  in  the  corresponding 
oil-flow  picture.  The  presence  of  two  'peaks'  may  correspond  to  the 
reattachment  and  separation  lines  caused  by  a  'separation  vortex'  which  was 
suggested  by  Token  with  his  'vortex  dominated  flow  field  model'  as  illustrated 
in  Fig. 10a. 

Very  recently.  Hung  et  al.  (Refs. 21  -  22)  reported  a  'numerical  solution  of 
three-dimensional  shock  wave  and  turbulent  boundary-layer  interaction' .  This 
solution  showed  good  agreement  with  some  experimental  data  at  the  Mach  number 
of  6,  as  shown  in  Fig. 11,  except  for  the  static  pressure  distribution  near 
the  shock  generator  when  the  shock  was  very  strong.  They  remarked  that 
further  study  is  needed  to  understand  the  details  of  viscous  interaction,  the 
problem  of  three-dimensional  separation  being  quite  different  from  the  two- 
dimensional  phenomenon.  They  also  pointed  out  that  three-dimensional 
turbulence  modelling  is  a  key  issue  for  further  numerical  simulation. 


1.2.  The  Present  Experimental  Study 

The  present  paper  attempts  to  build  a  flow  field  model  of  the  three- 
dimensional  glancing  shock/turbulent  boundary-layer  interaction  region,  on 
the  strength  of  (i)  oil-flow  patterns  (on  both  the  side  wall  and  the  shock 


generator  surface);  (ii)  vapour-screen  pictures;  (iii)  surface  static 


pressure  distributions  and  (iv)  heat  transfer  distributions.  The 
experimental  results  were  obtained  at  the  Cranfield  Institute  of  Technology 
with  two  separate  test  programmes  conducted  in  the  2.S  x  2.5  inch 
intermittent  tunnel'  and  the  9x9  inch  cdntinuoufe-tu'nhing’fuhhel,  'at  a 
Mach  number  of  approximately  2.5.  Some  of  the  experimental  data  from  the 
intermittent  tunnel  has  already  been  reported  by  Dickman  (Ref. 23), 

Davenport  (Ref. 24)  and  Kubota  (Ref. 25). 

The  new  flow  field  model  of  the  interaction  region  consists  of  two  different 
viscous  layers,  (1)  the  side-wall  boundary  layer  and  (2)  the  'induced 
layer'  originating  near  the  shock  generator  root  and  crossing  the  path  of 
the  side-wall  boundary  layer.  The  flow  field  model  leads  to  the  following 
hypothesis.  As  long  as  surface  stream  lines  of  the  side-wall  boundary 
layer  are  pliable  enough  to  be  'bent'  along  the  edge  of  the  induced  layer, 
no  separation  appears  there  even  if  the  surface-flow  deflection  exceeds  the 
shock  angle.  An  'ordinary  separation'  takes  place  when  the  induced  layer 
forces  the  surface  stream  lines  to  deflect  beyond  a  maximum  permissible 
angle  depending  on  the  side-wall  boundary-layer  properties.  The  surface 
stream  lines  then  lift  off  from  the  wall  surface,  sliding  over  the  induced 
layer  and  become  involved  in  a  vortical  free  layer,  as  shown  in  Fig. 10b. 

The  maximum  permissible  angle  may  be  determined  by  a  simple  angular  momentum 
relation. 

The  'ordinary  separation'  in  the  present  paper  is  defined  by  limiting  stream 
lines  leaving  the  surface  at  a  particular  line  (separation  line)  lying 
on  the  side  wall.  This  separation  will  be  marked  by  surface  stream  lines 
converging  and  coalescing  into  a  single  line,  which  can  be  experimentally 
detected  using  surface-flow  visualisation,  (oil-flow  pictures) . 
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2.  EXPERIMENTAL  PROGRAMMES 

Experimental  data  presented  in  this  paper  were  taken  in  two  separate  wind 
tunnel  test  programmes ,  at  almost  the  same  Mach  number  of  approximately  2.5. 

The  main  difference  between  these  two  test  programmes  was  in  the  scale 
of  the  shock  boundary- layer  interaction  region  to  be  investigated. 

A  small-scale  test  programme  was  conducted  in  the  2.5  x  2.5  inch  intermittent 
supersonic  wind  tunnel  (Fig. 12).  The  results  of  this  programme  have  been 
partially  reported  by  Dickman  (Ref. 23),  Davenport  (Ref. 24)  and  Kubota  (Ref. 25). 
The  large-scale  test  programme  was  undertaken  in  the  9x9  inch  continuous- 
running  supersonic  wind  tunnel  (Fig. 13)  and  some  data  obtained  in  it  have 
been  presented  by  Kubota  (Ref. 25). 

In  both  the  test  programmes,  the  interaction  region  was  provided  by  locating 
a  wedge  in  the  tunnel  working  section  so  that  an  oblique  shock  generated 
by  it  could  glance  across  the  turbulent  boundary  layer  growing  on  the 
tunnel  side  wall. 


2.1  2.5  x  2.5  inch  Wind  Tunnel  Test  Programme 

2.1.1.  Wind  Tunnel  Facility  and  the  Test  Conditions 

The  intermittent  wind  tunnel  has  a  cross  section  of  2.5  inches  by  2.5  inches 
at  its  working  section,  and  a  two-dimensional  unsysmetrical  supersonic 
nozzle  (Fig. 14).  By  interchanging  supersonic  nozzles  a  range  of  the  test 
Mach  number  from  1.6  to  4.3  is  available. 

The  wind  tunnel  is  ’driven'  by  air  at  atmospheric  pressure  sucked  down  into 
an  evacuated  tank.  The  low  pressure  tank  has  a  volume  of  approximately 
1600  cubic  feet  and  can  be  evacuated  by  a  25  kilowatt  electrically  driven 


pump. 
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Stagnation  conditions  for  this  tunnel  are,  therefore,  nearly  atmospheric 
and  tend  to  fluctuate  from  day  to  day  with  changing  weather. 


Stagnation  Pressure  and  Test  Mach  Number 

The  stagnation  pressure  was  measured  using  a  Pitot  probe  set  up  in  the 
tunnel  settling  chamber  and  connected  to  a  mercury  manometer.  This  pressure 
was  approximately  50  mmHg  lower  than  atmospheric  due  to  the  pressure  drop 
across  the  silica-gel  drying  system.  Throughout  this  test  programme,  the 
stagnation  pressure  was  710  ±  10  mmHg. 

The  static  pressure  was  measured  via  a  tapping  suitably  placed  in  the  top 
liner  of  the  tunnel. 

These  measured  stagnation  and  free-stream  static  pressures  determined  the 
free-stream  Mach  number,  through  the  isentropic  flow  relations.  In  this 
test  programue  which  used  a  'Mach  2.5'  nozzle  section,  the  test  Mach 
number  was  calculated  to  be  2.41  ±  0.03.  The  Schlieren  system  showed  that 
the  flow  was  steady  during  the  30  seconds  of  useful  running  time. 


Stagnation  Temperature 

The  stagnation  temperature  was  measured  by  means  of  a  Chromel-Almel  thermo¬ 
couple  located  in  the  tunnel  settling  chamber.  This  thermocouple  was 
connected  with  a  thermocouple  multimeter.  The  stagnation  temperature  for 
this  test  programme  was  in  a  range  of  290  ±  10  degrees  Kelvin. 

It  was  also  found  that  the  stagnation  temperature  tended  to  decrease  by  3  to 
4°K  during  the  tunnel  run.  This  temperature  behaviour  is  thought  to  be  due 


to  the  admittance  of  air  from  the  atmosphere,  through  a  large  volume  of 
inlet  piping  and  the  silica-gel  drying  system. 
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Reynolds  Number 

Isentropic  flow  relations  determined  the  free-stream  temperature  for  the 

measured  stagnation  temperature.  Southerland's  law  was  then  used  to 

calculate  the  free-stream  viscosity,  while  the  equation  of  state  provided 

the  air  density.  Thus,  after  measuring  the  free-stream  Mach  number,  the 

unit  free-stream  Reynolds  number  was  estimated  to  be  approximately  2.5  x  10~* 
per  inch. 


2.1.2  Test  Configuration  (See  Fig. 15) 

Shock  Generator 

To  generate  the  shock  to  interact  with  the  test  boundary  layer,  a  variable- 
angle  wedge  was  mounted  on  a  quadrant  suspended  from  the  tunnel  top  liner. 
The  pivot  point  of  the  wedge  was  designed  to  coincide  exactly  with  its 
leading  edge  resting  on  the  top  liner  surface.  This  configuration  made 
it  possible  to  increase  the  wedge  angle  up  to  approximately  10  degrees 
without  the  wind  tunnel  choking. 

The  wedge  had  a  chord  of  3.5  inches  and  spanned  the  tunnel  test  section. 
Sealing  grooves  were  cut  on  both  sides  of  the  wedge,  in  order  to  prevent 
pressure  leakage  between  the  wedge  and  the  side  walls. 

The  trailing  edge  of  the  wedge  generated  an  expansion  fan, and  the  oblique 
shock  from  the  wedge  leading  edge  reflected  from  the  tunnel  bottom  liner 
and  interacted  with  the  boundary  layer  there.  These  phenomena  limited  the 
effective  test  area.  The  flow  field  was  therefore  restricted  to  the  region 
upstream  of  the  wedge  trailing  edge  and  some  distance  away  from  the  tunnel 


-  11  - 


bottom  liner;  i.e.  a  zone  approximately  4  inches  in  length  and  2  inches 
in  width. 

On  the  top  surface  the  oblique  shock  did  not  occur  exactly  at  the  wedge 
leading  edge,  but  approximately  0.15  inches  upstream  of  it,  (according  to 
Schlieren  photographs)  presumably  due  to  an  influence  of  the  top-liner 
boundary  layer.  It  must  be  stressed  that  no  two-dimensional  flow  separation 
was  observed  at  the  wedge-top  liner  'corner'  with  the  angles  used  in  this 
experiment. 


Test  Boundary  Layer 

The  test  boundary  layer  was  a  fully  turbulent  boundary  layer  growing  along 
the  tunnel  side  wall.  The  boundary  layer  was  surveyed  by  using  Pitot  and 
static-pressure  probes  (described  in  the  next  chapter)  connected  to  a 
multi-tube  mercury  manometer.  These  surveys  were  made  at  five  stations  along 
a  swept  line  lying  2  inches  upstream  of  the  oblique  shock  at  the  shock 
generator  angle  of  5  degrees  (Fig.l6a).  The  survey  stations  were  selected 
so  as  to  give  the  undisturbed  boundary-layer  properties  in  almost  the  same 
distance  from  the  shock. 

The  Pitot-probe  surveys  showed  that  static  pressures  in  the  side-wall 
boundary  layer  were  almost  constant  (within  i  32)  and  equal  to  the  external- 
flow  static  pressure  and  the  free-stream  static  pressure  measured  on  the 
tunnel  top  liner. 

The  measured  static  and  total  pressure  distributions  in  the  boundary  layer 
enabled  the  density  and  velocity  profiles  to  be  calculated  by  the  method 
described  in  Appendix  A.  The  boundary- layer  properties  and  the  velocity 


profiles  are  shown  in  Figs. 16b  and  16c.  The  boundary  layer  thickness 
(defined  by  the  point  where  0.995  of  the  external-flow  velocity  is  reached) 
was  typically  0.2  inches.  The  corresponding  displacement  and  momentum 
thicknesses  were  estimated  to  be  approximately  0.07  inches  and  0.015  inches 
respectively. 


Modification  of  the  Test  Boundary  Layer 

For  some  of  the  tests  an  attempt  was  made  to  increase  the  side-wall  boundary- 
layer  thickness  by  using  vortex  generators,  roughness  elements  and  air 
injection.  Many  configurations  were  tried  at  various  locations  in  the  tunnel. 
The  most  effective  type  of  disturbance  was  a  group  of  small  poles  (0.5 
inches  in  height,  0.5  inches  in  length  and  0.1  inches  in  width)  stuck  on 
the  tunnel  side  wall  a  little  upstream  of  the  supersonic  nozzle  throat 
(Fig. 17).  This  type  of  roughness  increased  the  boundary  layer  thickness 
to  approximately  0.3  inches.  Meanwhile,  the  corresponding  displacement 
and  momentum  thicknesses  were  increased  to  0.085  and  0.018  inches 
respectively.  These  modified  boundary-layer  properties  are  plotted  in  Fig. 16b 
together  with  the  original  properties. 

This  device  also  'flattened'  the  boundary-layer  thickness  distribution 
over  the  flow  field,  as  shown  in  Fig. 16b.  It  seems  that  this  device 
invigorated  the  side  wall  boundary  layer  and  made  it  more  resistant  to 
change . 


2.1.3  Test  Techniques 


With  the  test  configuration  using  the  original  side-wall  boundary  layer 
(without  the  roughness  elements)  the  following  test  techniques  were  applied: 

(1)  Schlieren,  normal  to  the  tunnel  side  wall 

(2)  Oil-flow  on  the  side  wall  and  on  the  shock  generator  surface. 

(3)  Vapour-screen  in  cross  sections  normal  to  the  free  stream. 

(4)  Static  pressure  distribution  measurements  on  the  side  wall. 

(5)  Heat  transfer  measurements  on  the  side  wall,  by  using  a  quasi¬ 

transient  'slug  calorimeter'  technique. 

(6)  Encapsulated  liquid  crystal  coatings  to  give  surface  temperature 

patterns  on  the  side  wall. 

(7)  Viscous  layer  surveys  using  Pitot  probes  and  a  yaw  meter. 

The  results  obtained  with  these  techniques  in  this  test  configuration  will 
be  specified  as  pertaining  to  Model  A. 

In  the  modified  test  boundary  layer,  the  following  techniques  were  applied: 

(1)  Schlieren 

(2)  Oil-flow 

(3)  Vapour-screen 


The  results  will  be  specified  as  referring  to  Model  AM 
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2.2  9x9  inch.  Wind  Tunnel  Test  Programme 

2.2.1  Wind  Tunnel  Facility  and  the  Test  Conditions 

The  continuous-running  tunnel  has  a  cross  section  of  9  inches  by  9  inches 
at  the  working  section,  and  a  two-dimensional  supersonic  nozzle  (Fig. 18). 

The  supersonic  nozzle  is  also  interchangeable  to  cover  a  range  of  Mach 
number  from  1.7  to  2.5.  This  tunnel  is  driven  by  a  centrifugal  compressor 
with  two  373  kilowatt  electric  motors. 

The  stagnation  pressure  is  automatically  maintained  at  a  pre-set  value 
in  the  range  from  0.7  to  14.7  p.s.i.a.  The  stagnation  temperature  could 
be  controlled  in  the  range  approximately  from  5°K  lower  to  10°K  higher  than 
the  ambient  temperature,  by  regulating  the  amount  of  water  in  the  tunnel 
cooling  system.  However,  it  is  difficult  to  control  the  temperature 
accurately  because  the  adjustment  involves  a  long  time  constant  (approximately 
from  1.5  to  2.5°K/Min.). 


Stagnation  Pressure  and  Test  Mach  Number 

For  this  test  programme,  the  stagnation  pressure  was  preset  to  be 
approximately  4  p.s.i.a.  The  stagnation  pressure  was  monitored  from  the 
pressure  gauge  mounted  in  the  tunnel  control  panel.  When  higher  accuracy 
of  the  measured  stagnation  pressure  was  required,  such  as  in  the  boundary 
layer  surveys,  a  Pitot  probe  set  up  in  the  tunnel  settling  chamber  was 
used  to  provide  it  through  the  diaphragm- type  pressure  transducer  of  a 
scanning  valve  system  (described  in  the  next  chapter) .  The  stagnation 
pressure  was  estimated  to  be  4.0  ±  0.2  p.s.i.a. 
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The  free-stream  static  pressure  was  obtained  on  the  tunnel  top-liner 
surface  at  two  stations  located  between  the  supersonic  nozzle  throat  and 
the  test  region  (Fig. 18).  The  pressures  measured  on  the  top  liner  were 
confirmed  to  be  almost  equal  to  the  real  free-stream  static  pressure  during 
the  survey  through  the  boundary  layer. 

These  measured  stagnation  and  free-stream  static  pressures  gave  the  test 
Mach  number,  in  the  same  way  as  in  the  2.5  x  2.5  inch  wind  tunnel  test 
programme.  The  test  Mach  number  was  calculated  to  be  2.30  ±  0.05. 


Stagnation  Temperature 

To  measure  the  stagnation  temperature,  a  Chrome 1-Almel  thermocouple  probe 
was  mounted  in  the  tunnel  settling  chamber  and  continuously  indicated  the 
stagnation  temperature  during  the  tunnel  running  time. 

For  the  heat  transfer  measurements,  it  was  required  to  maintain  a  low 
stagnation  temperature  to  provide  a  large  difference  between  the  wall 
temperature  and  the  air  recovery  temperature.  In  this  case,  the  stagnation 
temperature  was  kept  between  282  and  285  degrees  Kelvin,  by  fully  operating 
the  tunnel  cooling  system.  For  other  measurements,  the  temperature  lay  in 
the  range  from  283  to  293  degrees  Kelvin. 


Reynolds  Number 

In  the  same  way  as  in  the  2.5  x  2.5  inch  wind  tunnel  test  programme,  the 

5 

unit  free-stream  Reynolds  number  was  estimated  to  be  approximately  0.8  x  10 
per  inch,  for  the  stagnation  pressure  of  4  p.s.i.a. 
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2.1.2.  Test  Configuration  (See  Fig. 19) 

Shock  Generator 

A  thin  wedge  was  mounted  on  a  rotating  disk  (9  inches  in  diameter)  mounted 
in  the  tunnel  side  wall.  This  wedge  was  made  from  a  J"flat  plate  with  the 
leading  edge  chamfered  to  an  angle  of  15°.  The  chord  was  6"  and  the  span  7" 
so  that  the  model  occupied  about  80Z  of  the  working  section  width. 

The  shock  generator  angle  was  pre-set  by  turning  the  disk.  The  pivot  point 
did  not  coincide  with  the  wedge  leading  edge  hence  the  wedge  leading  edge 
moved  with  changing  shock  generator  angle.  The  movement  was  approximately 
0.8  inches  in  both  the  directions  parallel  and  normal  to  the  tunnel  liner, 
for  an  angle  change  of  0  to  15  degrees.  It  was  checked  that  the  space 
between  the  wedge  leading  edge  and  the  tunnel  bottom  liner  was  always  larger 
than  1  inch  in  order  to  avoid  interference  from  the  boundary  layer  growing 
along  the  tunnel  bottom  liner. 

This  experimental  configuration  provided  a  test  flow  area  of  approximately 
6  inches  by  6  inches,  without  any  interference  from  shock  reflections, 
trailing  edge  and  the  tunnel- liner  boundary  layers. 


Test  Boundary  Layer 

The  test  boundary  layer  was  the  fully  developed  turbulent  boundary  layer 
growing  along  the  tunnel  side  wall.  This  boundary  layer  was  surveyed  using 
pitot  and  static-pressure  tubes  (with  the  same  dimensions  as  in  the  previous 
test  programme)  connected  to  the  pressure  transducer  of  the  scanning  valve 
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system.  Boundary  layer  surveys  were  made  without  the  wedge  at  three  stations. 
These  stations  were  located  along  a  line  normal  to  the  free  stream,  4.1 
inches  upstream  of  the  rotatory-disk  centre  (Fig. 20a). 

The  survey  showed  that  the  static  pressure  in  the  side-wall  boundary 
layer  was  constant  and  equal  to  the  external-flow  static  pressure  and  the 
free-stream  static  pressure  measured  on  the  tunnel  top  liner,  (to  within  ±  5Z) . 

The  boundary-layer  properties  and  the  velocity  profiles  are  presented  in 
Figs. 20b  and  20c.  The  boundary  layer  thickness  (defined  by  a  point  where 
0.995  of  the  external-flow  velocity  is  reached)  was  approximately  0.62  inches, 
and  the  corresponding  displacement  and  momentum  thicknesses  were  estimated 
to  be  0.15  and  0.040  inches  respectively.  These  results  confirm  that  the 
side-wall  boundary  layer  was  almost  uniform  over  the  test  area. 


2.2.3.  Test  Techniques 

In  this  test  programme ,  the  following  test  techniques  were  applied: 

(1)  Schlieren,  normal  to  the  tunnel  side  wall. 

(2)  Oil-flow  visualization  on  the  side  wall  and  on  the  shock 

generator  surface. 

(3)  Vapour-screen  pictures  of  cross  sections  normal  to  the  free  stream. 

(4)  Static  pressure  distribution  measurements  on  the  side  wall. 

(5)  Heat  transfer  measurements  on  the  side  wall,  using  a 

quasi-transient  'thin- film-gauge  calorimeter*  technique. 


The  results  obtained  with  these  techniques  in  this  test  programme  will  be 
specified  as  referring  to  Model  B. 
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3.  EXPERIMENTAL  TECHNIQUES  AND  THE  ARRANGEMENTS 
3.1.  Schlieren  Techniques 

Schlieren  pictures  were  taken  for  all  model  configurations  to  check  that 
the  test  flow  pattern  was  as  required.  The  pictures  were  taken  normal  to 
the  tunnel  side  walls  and  the  knife-edge  was  adjusted  so  as  to  give  good 
contrast. 

For  Models  A  and  AM,  the  configuration  of  the  shock  generator  (suspended 
from  the  tunnel  top  liner)  allowed  glass  windows  to  be  fitted  on  both  sides 
of  the  tunnel  test  section.  Hence  a  conventional  single  pass  Schlieren 
system  was  used  and  this  gave  good  quality  pictures. 

In  the  case  of  Model  B  the  wedge  was  mounted  on  an  aluminium  disk,  free  to 
rotate  within  the  tunnel  side  wall,  made  from  steel.  Hence  a  glass  window 
could  be  fitted  on  one  side  only.  Therefore  the  disk  surface  was  mechanically 
polished  to  be  as  smooth  as  possible  and  a  double-pass  Schlieren  system 
was  used.  Although  this  optical  system  could  not  escape  some  interference 
due  to  the  roughness  of  the  disk  surface,  the  Schlieren  pictures  were  adequate. 


3.2.  Oil-Flow  Pictures 

To  visualize  the  surface  flow  on  the  tunnel  side  wall  and  on  the  shock 
generator  surface  (Fig. 21),  a  mixture  of  titanium  dioxide  suspended  in  motor 
oil  was  used.  A  drop  of  oleic  acid  was  added  to  the  mixture  to  prevent 
coagulation  of  the  titanium  dioxide. 

The  surfaces  to  be  painted  with  the  white  oil  mixture  were  prepared  with  a 
matt-black  finish  to  provide  the  greatest  photographic  contrast. 
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To  avoid  smeared  oil  traces  of  the  surface  flow  patterns  (caused  by 
stopping  the  tunnel)  the  pictures  were  taken  through  the  glass  windows  while 
the  tunnel  was  running,  except  for  some  pictures  of  the  shock-generator 
surface  flow  patterns. 


3.3  Vapour-Screen  Pictures 

To  visualize  the  flow  features  in  cross  sections  of  the  test  flow  field 
normal  to  the  free  stream  (Fig. 22),  the  vapour-screen  technique  was  applied. 
The  principle  of  this  technique  is  as  follows. 

As  the  air  expands  through  the  supersonic  nozzle  into  the  tunnel  working 
section  and  cools,  the  moisture  in  the  air  tends  to  condense  out  forming  a 
fog.  The  fog  particles  can  scatter  light  in  proportion  to  the  number  of 
particles  per  unit  volume.  Therefore,  when  illuminated  by  a  narrow  'sheet' 
of  light  they  build  up  a  'vapour- screen'  which  indicates  disturbed  fog 
caused  from  disturbances  in  the  flow  field.  (Details  are  given  in  the  paper 
by  McGregor ,  Ref . 26) . 

In  the  application  of  this  technique  to  the  present  experiments,  two 
important  questions  arose.  They  were  how  to  control  the  fog  density  and  how 
to  increase  the  contrast  of  the  vapour-screen  pictures. 

For  Models  A  and  AM  tested  in  the  intermittent  tunnel,  a  small  steam-boiler 
with  an  electric  heater  was  used  to  supply  enough  fog  in  the  working  section. 
The  steam  was  injected  into  the  tunnel  settling  chamber.  After  filling 
the  settling  chamber  with  the  steam,  the  tunnel  was  started.  This  procedure 
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provided  enough  time  in  which  to  take  the  pictures  (approximately  5  seconds). 
However,  condensation  of  the  steam  also  formed  water  droplets  on  the  cooled 
window  glass  and  tended  to  spoil  the  pictures  taken.  The  glass  window 
had  to  be  cleaned  after  every  run. 

In  the  case  of  Model  8,  the  continuous-running  tunnel  made  it  easy  to 
control  the  density  of  the  fog.  By  injecting  a  given  amount  of  water  into 
the  tunnel  settling  chamber,  a  sufficient  density  of  the  fog  was  obtained 
for  long  enough  to  take  a  picture. 

In  practice  it  was  difficult  to  determine  quantitatively  the  amount  of  injected 
steam  or  water  to  give  the  best  vapour-screen  pictures.  The  fog  density 
depended  critically  on  the  ambient  conditions.  By  eye  it  was  was  hard  to 
distinguish  the  blurred  vapour-screen  images.  Therefore,  in  most  cases 
contrast  of  the  photographs  taken  needed  to  be  improved.  This  improvement 
was  done  photographically  by  recopying  or  reprinting  the  pictures.  However, 
this  photographic  process  tended  to  emphasize  any  marks  or  scratches  on 
the  pictures. 

In  the  present  programmes,  the  illumination  used  was  a  250  W  slide-projector. 

A  loaded  slide  with  a  small  slit  produced  a  sheet  of  light  which  was  further 
narrowed  using  a  convex  lens.  The  sheet  was  approximately  0.2  inches  thick 
across  the  test  section.  The  pictures  were  taken  on  Kodak  Trix-pan  film 
(ASA  400)  with  an  exposure  time  of  1  to  4  seconds  for  an  aperture  of  f/4. 

The  film  developer  was  Paterson  ACUSPEED  FX-20  which  can  increase  the 
effective  film  speed  by  up  to  three  times.  The  recopy  and  reprint  of  the 
pictures  employed  ILFORD  PANF  film  (ASA  50)  and  ILFOSPEED  4.1  M  printing 
paper  respectively. 
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3.4.  Static  Pressure  Distribution  Measurements 

The  static  pressure  distributions  were  measured  by  means  of  holes  (0.02 
inches  in  diameter)  arranged  over  the  side  wall  of  the  2.5  x  2.5  inch  tunnel 
and  on  the  metal  disk  of  the  9x9  inch-tunnel  side  wall.  The  distributions 
of  these  holes  in  the  two  tunnels  are  shown  in  Fig. 23a  and  Fig. 23b. 

On  the  side-wall  of  the  2.5  x  2.5  inch  tunnel,  two  hundred  and  ten  pressure 
tappings  were  arranged  in  seven  rows  of  30  tappings.  On  the  metal  disk 
of  the  9x9  inch  tunnel  wall  one  hundred  and  fifty-six  orifices  were 
arranged  in  6  rows. 

The  measured  pressures  were  indicated  by  a  digital  voltmeter  on  the  control 
panel  of  the  scanning  valve  system  feeding  a  diaphragm- type  pressure 
transducer  (Fig. 23c).  The  system  was  calibrated  with  the  known  ambient 
pressure  and  a  low  pressure  (approximately  0.05  mmHg)  provided  by  a  vacuum 
pump,  for  each  individual  scanning  port  (30  stations  for  Model  A  and  93 
stations  for  Model  B) . 

Particularly  for  Model  A,  tested  in  the  intermittent  tunnel,  the  tunnel 
running  time  of  less  than  30  seconds  was  too  short  to  measure  the  correct 
pressures.  This  was  because  the  lines  connecting  the  tappings  to  the 
scanning  valve  system  had  a  large  volume  which  filled  too  slowly  to  allow 
the  measured  pressures  from  reaching  their  final  values.  This  pressure 
delay  was  improved  by  using  a  guillotine  clamp  to  enable  accumulation 
of  the  pressures  through  repeated  tunnel-running.  The  clamp  allowed  the 
pressure  to  build  up  in  the  lines  when  the  tunnel  was  running  and  trap  them 
when  the  tunnel  stopped.  This  procedure  was  repeated  several  times  until 
the  measured  pressures  reached  their  equilibrium  state. 
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3.5.  Heat  Transfer  Measurements 

To  measure  the  heat  transfer  distributions  a  quasi-transient  calorimeter 
technique  was  applied.  Both  the  experimental  procedure  and  the  type  of 
calorimeter  employed,  were  different  between  the  two  tunnel  test  programmes. 
This  was  because  of  a  basic  difference  between  the  two  sets  of  test 
conditions,  i.e.  intermittent  and  continuous. 

In  the  2.5  x  2.5  inch  tunnel  test  programme,  a  modified  Princeton  University 
slug-calorimeter  technique  was  used.  Meanwhile,  for  the  9x9  inch  tunnel 
test  programme,  a  thin-f ilm-gauge  calorimeter  technique  was  devised  and 
developed. 

The  principles  of  these  two  techniques  are  discussed  in  Appendix  B. 

3.5.1.  A  Modified  Princeton  University  Slug  Calorimeter  Technique 

The  original  technique  using  the  Princeton  University  type  of  slug 
calorimeter  is  fully  described  by  Oskam  et  al.  in  Ref. 27.  The  modified 
technique  is  detailed  by  Davenport  in  Ref. 24.  Some  discussion  of  these 
techniques  is  given  in  Appendix  B2. 

The  Original  Technique 

When  a  piece  of  metal  (slug)  of  known  characteristics  is  exposed  to  a 
thermal  medium,  the  temperature  variation  of  the  slug  with  time  can  give 
the  heat  transfer  rate  between  the  slug  surface  and  the  medium  directly. 

The  Princeton  University  type  of  calorimeter  employs  a  copper  slug  which 
can  be  heated  by  a  tiny  jet  of  hot  air  and  insulated  from  the  metal  container 
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by  a  nylon  phenolic  material.  When  the  tunnel  is  started  the  tiny  jet  of 
hot  air  is  supplied  to  maintain  the  slug  at  the  temperature  of  the  wall 
surrounding  it,  which  is  higher  than  the  adiabatic  wall  temperature  (or 
recovery  temperature)  of  the  test  flow.  When  the  data  acquisition  system 
is  ready,  extra  heat  is  added  to  the  slugs  to  raise  their  temperature , 
typically  by  20°K  (see  Ref. 27).  When  the  heating  jet  is  stopped,  the  data 
acquisition  system  is  started  and  the  slugs  are  allowed  to  cool  naturally, 
down  through  the  wall  temperature  to  eventually  reach  the  adiabatic  wall 
temperature,  (Fig. 24a).  By  evaluating  the  gradient  of  the  slug-temperature 
time  history  at  the  point  where  the  slug  temperature  crosses  the  wall 
temperature,  the  heat  transfer  rate  from  the  wall  to  the  air  can  be 
determined  (see  Appendix  B2) . 


The  Modified  Technique 

In  the  modified  technique,  the  "tiny  jet  of  hot  air"  was  not  used.  Instead 
the  temperature  difference  to  drive  the  slug  calorimeter  was  provided  by 
the  difference  between  the  wall  temperature  and  the  adiabatic  wall 
temperature.  This  was  done  because  it  was  found  that  heat  convection  during 
the  heating  process  disturbed  the  test  wall  surface  and  hence  affected 
neighbouring  gauges. 

The  experimental  procedure  was  as  follows.  The  slug  calorimeters,  as  well 
as  the  tunnel  side  wall,  were  left  to  reach  room  temperature.  When  the 
intermittent  tunnel  was  started,  the  slugs  were  cooled  by  the  air  flow  down 
to  the  adiabatic  wall  temperature.  The  temperature-time  history  is  shown 
in  Fig. 24b  The  slugs  have  much  smaller  thermal  capacities  than  the  tunnel 
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wall  and  hence  they  cooled  more  quickly.  Therefore,  the  only  time  that 
the  8 lug  temperatures  coincided  with  the  wall  temperature  was  at  t  •  0; 
i.e.  at  the  moment  that  the  tunnel  was  started.  Hence  the  heat  transfer 
rates  were  derived  from  the  temperature-time  plots  by  measuring  the  slope 
at  t  ■  0. 


Test  Surface  Installation 

To  accommodate  the  slug  calorimeters  in  the  tunnel  side  wall,  an  aluminium 
plug  of  2  inches  diameter  was  used,  having  three  insulated  slugs  on  its 
flat  surface  (Fig. 25a).  Two  wall  temperature  sensors  (Chromel-Alumel 
thermocouples  embedded  in  copper)  were  also  flush  with  the  surface.  The 
plug  was  kindly  loaned  by  Princeton  University  and  had  been  used  in  an 
earlier  study  by  Oskam  et  al.  (Ref. 27). 

To  mount  the  plug  in  the  test  section,  an  eccentric  rotatable  system  was 
fabricated.  A  4  inch  diameter  plug  was  fitted  into  a  hole  of  equal  diameter 
in  the  tunnel  side  wall.  In  the  4  inch  diameter  plug,  a  2  inch  diameter  hole 
was  drilled  C.75  inches  off  centre.  The  2  inch  diameter  plug  holding  the 
slug  calorimeters  was  fitted  in  this  hole.  The  whole  system  is  illustrated 
in  Fig. 25b.  Rotation  enabled  a  wide  range  of  measuring  stations  to  be  covered. 


Data  Acquisition  System  and  Data  Reduction 

Chromel-Alumel  thermocouples  embedded  in  the  slugs  and  the  plug  container 
gave  the  temperature  data  to  be  used  in  calculating  the  heat  transfer  rates. 
The  thermocouple  leads  were  joined  to  a  series  of  linear  amplifiers  and  the 
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amplified  signals  were  recorded  by  an  ultra-violet  recorder.  This  measurement 
system  was  calibrated  using  the  water-bath  method,  in  advance  of  the  tunnel 
test. 


When  the  tunnel  was  running,  the  slug  temperatures  dropped,  approaching  the 
adiabatic  wall  temperature  exponentially.  In  contrast  the  tunnel  wall 
temperature  changed  very  little  because  of  its  large  thermal  capacity. 
However,  it  was  also  observed  that  the  freestream  stagnation  temperature 
decreased  by  approximately  3  to  4°K,  in  the  first  12  seconds  of  the  running 
time.  This  temperature  behaviour  was  because  initially  the  tunnel  sucks 
through  the  air  that  has  been  filling  the  intake  duct  and  settling  chamber. 
Subsequently  the  tunnel  draws  in  the  colder  fresh  air  from  outside  the 
laboratory.  Because  of  this  temperature  variation  the  measured  slug 
temperatures  had  to  be  corrected,  as  shown  in  Fig. 26. 


From  the  corrected  temperature  time-history,  the  temperature  gradient  at 
t  -  0  was  evaluated.  This  gradient  gave  the  heat  transfer  rate,  q,  from 
the  relation: 


q 


(3  -  1) 


where  Tg  is  the  slug  temperature  while  Cg  and  Ag  are  respectively 
thermal  capacity  and  surface  area  of  the  slug  calorimeter. 


The  adiabatic  wall  temperature,  T^  ,  is  calculated  from  the  measured 
stagnation  temperature,  Tq  ,  and  the  measured  free  stream  Mach  number, 
Mj  by 
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Ta„/To  ’  (1  +*Ii!mi2)/(1  +ITmi2)  >  O  -  2) 

where  r  end  y  ere  respectively  the  recovery  factor  and  the  ratio 
of  specific  heats,  which  may  be  taken  as  0.89  (for  a  turbulent  boundary 
layer)  and  1.4.  In  the  present  test  conditions,  the  ratio  of  the  wall 
temperature,  T  ,  to  the  adiabatic  wall  temperature  can  be  written  as 

T  /T  -  1.06  T  /T  (3  -  3) 

w  aw  w  o  v 

and  the  temperature  difference  to  drive  the  slug  calorimeter  is 
approximately 

T  -  T  -  18°K  .  (3-4) 

w  aw  v 

Using  the  slug  calorimeter  technique,  the  heat  transfer  to  the  side- 
wall  with  no  shock  generator  (bare  plate  heat  transfer)  was  measured  as 
shown  in  Fig. 27.  In  the  figure,  the  local  Stanton  number  is  defined  by 

ch  *  i*!  ui  CP  (T.  -  T.„>  «  -  5> 

where  Cp  is  specific  heat  of  the  air  while  and  u^  are  respectively 
the  density  and  the  velocity  of  the  free  stream.  The  results  show  values 
approximately  10Z  higher  than  a  prediction  based  on  Van  Driest  II  theory 
(Ref. 31). 

3.5.2  Thin-Film-Gauge  Calorimeter  Technique 

Details  of  the  thin-filar  gauge  calorimeter  technique  are  described  in 
Appendix  B3,  relating  to  the  quasi-transient  method  used  in  the  continuous- 
running  tunnel. 


Outline  of  the  Technique 


The  thin  film  gauge  is  usually  employed  as  a  thermometer,  especially  in 
short  duration  tunnels  with  a  high  ratio  of  stagnation  temperature  to  wall 
temperature  (T  /T  ) .  The  resistance  of  the  gauge  changes  almost  in  proportion 
to  its  temperature.  This  resistance  variation  is  calculated  by  operating 
the  gauges  at  a  constant  current  and  measuring  the  changing  voltage  drop 
across  the  film. 

If  the  intensity  of  the  input  electric  current  is  increased  then  the  gauge 
can  become  a  heater  and  locally  raise  the  temperature  of  the  substratum. 
Therefore  if  the  thermal  characteristics  of  the  gauge  and  its  substratum 
are  known,  the  gauge-substratum  combination  is  usable  as  a  calorimeter  with 
a  'built -in- heater  and  thermometer'. 

The  advantages  and  disadvantages  of  this  technique  as  compared  with  the 
previous  slug-calorimeter  technique  are  as  follows: 

(1)  Since  the  gauge  film  thickness  is  of  the  order  of  O.lym,  the 
thermal  resistance  between  the  'thermometer'  and  the  'surface'  is  extremely 
small.  This  makes  it  easy  to  obtain  a  very  high  sensitivity  calorimeter 
which  can  respond  to  small  changes  in  heat  transfer  and  which  has  a 

quick  response  time. 

(2)  The  thermal  capacity  of  the  gauge  is  much  smaller  than  that  of  the 
substratum.  If  the  substratum  is  used  as  the  test  surface,  the  temperature 
of  the  calorimeter  will  be  virtually  identical  to  that  of  the  test  surface. 

As  a  result,  the  measured  heat  transfer  is  expected  to  be  very  close  to 
the  real  heat  transfer  from  the  test  surface. 
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(3)  To  read  the  temperature  some  electric  current  must  pass  through  the 
’thermometer'.  This  electric  current  could  cause  some  heating  which 
disturbs  the  surface  temperature  to  be  measured.  To  reduce  this 
disturbance  to  a  negligible  level,  the  intensity  of  the  input  electric 
current  must  be  carefully  chosen,  especially  for  tunnels  with  a  low 
temperature  difference  (T  -T . 


Quasi-Transient  Technique  in  the  Continuous-Running  Tunnel 

The  quasi-transient  technique  was  used  for  the  measurements  made  in  the 
continuous  wind-tunnel.  The  tunnel  is  in  the  laboratory  and  the  walls 
are  normally  at  room  temperature.  The  tunnel  is  equipped  with  a  heat 
exchanger  (water  cooler)  to  keep  the  air  stagnation  temperature  constant 
during  the  run.  By  operating  the  cooler  at  maximum  conditions  it  is 
possible  to  hold  Tq  below  T^  for  a  considerable  time  because  of  the  large 
thermal  capacity  of  the  tunnel  test  section  walls.  The  tunnel  conditions 
were  monitored  by  the  stagnation  temperature ,  Tq  ,  and  the  reference  wall 
temperature  measured  on  the  metal  disk  surface,  .  Strictly  speaking  the 

stagnation  and  wall  temperatures  were  changing  slightly  with  time,  through¬ 
out  the  tunnel  run,  as  shown  in  Fig. 28.  In  the  figure,  the  temperatures  are 
plotted  against  the  gauge  station  numbers  measured  in  turn,  with  the 
interval  time  of  approximately  2  minutes  between  readings.  However,  as 
indicated  in  the  same  figure,  the  difference  between  the  reference  wall 
temperature  and  the  adiabatic  wall  temperature  (calculated  in  the  same  manner 
as  in  3.5.1.)  could  be  maintained  at  almost  a  constant  value  in  the  range 
from  16  to  18°K. 
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The  procedure  to  obtain  the  temperature  time  history  (Fig. 29)  is  as  follows. 
The  gauge  to  be  operated  was  heated  to  a  temperature  from  3  to  5°K 
higher  than  the  reference  wall  temperature.  This  was  done  by  passing  a 
steady  current  of  around  30mA  through  the  gauge.  When  the  gauge  temperature 
had  nearly  reached  its  new  temperature,  the  gauge  was  switched  over  to 
a  much  smaller  electric  current  of  3mA  to  measure  the  gauge  temperature. 

As  the  gauge- substratum  combination  (calorimeter)  was  cooled  by  the  air- 
stream,  the  indicated  gauge  temperature  exponentially  approached  the  local 
wall  temperature.  This  process  provided  the  temperature-time  history 
which  when  compared  with  the  corresponding  history  for  wind-off  conditions 
can  be  reduced  to  a  heat  transfer  coefficient  by  the  method  described  in 
Appendix  B3. 


Thin-Film-Gauge  Calorimeter  Installation 

The  calorimeter  itself  is  no  different  from  the  standard  thin  film  gauge. 

The  arrangement  of  the  63  gauges  used  is  shown  in  Fig. 30a. 

The  gauges  were  platinum  films  sintered  on  a  Pyrex  glass  plate  at  its 
surface.  (The  sintering  process  is  described  in  Refs.  28-30).  The  size 
of  these  gauges  was  approximately  0.2  inches  long  and  0.04  inches  wide. 

The  Pyrex  glass  plate  was  6  inches  long  and  5  inches  wide,  while  its  thickness 
was  0.058  ±  0.001  inches.  Through  this  glass  plate,  0.04  inch  diameter 
holes  were  drilled  at  each  end  of  all  the  gauges.  A  silver  paint  was  used 
to  draw  lead  lines  through  the  holes,  from  the  gauge  ends  to  the  back 
surface  of  the  glass  plate.  On  the  silver-paint  lines  of  the  back  surface, 
copper  leads  were  soldered.  The  Pyrex  plate  was  flush  mounted  in  the 
rotatable  metal  disk  of  the  tunnel  side  wall,  by  using  peisin  wax. 
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The  copper  leads  were  drawn  out  of  the  tunnel  aide  wall  through  a  small 
hole  on  the  disk. 

The  gauges  were  calibrated  with  a  silicon-oil  bath  method,  in  advance  of 
installing  in  the  metal  disk.  Although  the  resistance  of  the  63  gauges 
ranged  from  16  to  71  ohms,  the  temperature  coefficient  of  resistance  was 
measured  as  always  lying  in  the  range  2.13  ±  0.07 /°K  for  a  1.00mA  input 
electric  current. 

To  measure  the  reference  wall  temperature,  a  Copper-Constantan  thermocouple 
was  flush  mounted  on  the  metal  disk  surface,  upstream  of  the  gauge  plate. 


Installation  of  the  Electronics 

For  the  quasi-transient  calorimeter  technique,  an  electronic  system  was 
prepared.  The  block  diagram  is  shown  in  Fig. 30b,  together  with  data 
acquisition  system,  and  the  detailed  circuits  are  presented  in  Appendix  B5. 

The  electronic  system  consisted  of  three  principal  units,  the  Wheatstone- 
bridge  balance  circuit,  the  current  regulator  and  the  amplifier.  The 
Wheatstone  bridge  network  incorporated  the  usual  balancing  arms  so  that 
steady  signals  were  'removed'.  Hence  the  amplifier  only  receives  the  voltage 
differences  emanating  from  variations  in  gauge  resistance.  The  Wheatstone 
bridge  circuit  used  'thick  film'  resistors  to  obtain  high  temperature 
stability  and  a  low  level  of  noise.  The  current  regulator  was  designed  to 
supply  5  different  intensities  of  electric  current  in  the  range  from  0.3  to 
50.0mA,  by  handling  toggle  switches  on  the  control  panel  of  the  system. 
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This  was  necessary  because  the  gauges  manufactured  had  widely  differing 
resistances  which  therefore  required  different  electric  currents  for 
heating  the  calorimeters  and  measuring  the  local  surface  temperatures. 

The  regulated  current  varied  less  than  0.1Z  for  resistances  ranging  from 
10  to  150  ohms.  The  amplifier  was  built  up  from  two  operational-amplifier 
integrated  circuits.  The  minimum  gain  margin  per  amplifier  was  designed 
to  be  more  than  60  db,  to  provide  high  stability.  The  complete  amplifier 
circuit  gave  more  than  50  db  maximum  amplification  with  high  stability. 

In  all  of  this  electronic  system,  the  noise  level  converted  to  input  signal 
was  less  than  0.05yV,  and  the  zero-point  drift  converted  was  from  10  to 
2 Op V/ hour  for  normal  operation  after  30  minutes  warming  up. 


Data  Acquisition  System  and  Data  Reduction 

The  thin  film  gauges  gave  the  variation  of  local  wall-surface  temperatures 
with  time.  The  thermocouple  embedded  in  the  tunnel  side-wall  surface 
indicated  the  reference  wall  temperature.  The  tunnel  thermometer  system  with 
thermocouple  transducer  provided  the  stagnation  temperature.  The  signals 
of  these  temperatures  were  properly  amplified  with  linear  amplifiers,  and 
filtered  through  a  C-R  filter  with  a  time  constant  of  60  milliseconds. 

They  were  recorded  by  an  ultra-violet  recorder,  whilst  being  continually 
monitored  (by  means  of  a  digital  voltmeter  for  the  local  wall-surface 
temperatures,  a  thermocouple  multimeter  for  the  reference  wall  temperature 
and  the  tunnel  thermometer  system  for  the  stagnation  temperature) . 

When  the  calorimeter  of  the  gauge-substrate  combination  was  operated 
through  the  heating  process,  the  gauge  temperature  approaches  the  local  wall 
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surface  temperature  exponentially.  Some  typical  time  histories  of  the 
recorded  gauge-temperatures  are  shown  in  Fig. 29b.  As  discussed  in  Appendix  B3 
this  gauge-temperature,  T  ,  ought  to  behave  in  the  form  given  by 
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where  and  T  ^  are  the  local  wall-surface  temperature  and  initial 

gauge  temperature,  and  h^  and  h^  are  the  heat  transfer  coefficients 

to  the  air  and  to  the  wall  (heat  loss) .  Typical  measured  gauge-temperature 

data  treated  by  using  equation  (3-6)'  are  plotted  in  Fig. 31,  where 

T  .*  *  AT  .  ,  A  being  an  arbitrary  constant, 
gi  gi 

According  to  (3  -  6)',  the  temperature  data  ought  to  lie  on  a  straight  line 
(Fig. 31).  However,  the  measured  temperature  data  actually  deviate  from 
the  expectation,  in  the  first  2-5  seconds  after  switching  the  calorimeter 
system  from  heating  to  cooling.  This  temperature  behaviour  is  considered 
to  be  due  to  a  transitional  temperature  re-distribution  in  the  calorimeter 
(Appendix  B4) .  After  this  short  transition  period,  the  temperature  data 
followed  the  line  predicted  by  (3  -  6) ' . 

In  the  present  data  reduction,  the  following  formula  derived  from  (3  -  6) ’ 
was  employed: 
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The  sampling  interval  time.  At,  was  chosen  to  be  one  second. 

The  heat  loss  (h ^  was  measured  when  the  tunnel  was  not  running  (wind-off) . 
It  was  assumed  that  this  h^  was  not  changed  by  running  the  tunnel  (wind-on) . 


To  examine  the  influence  of  the  heating  process  on  the  temperature  data 
obtained,  some  measurements  in  the  wind-off  condition  were  conducted  while 
changing  the  intensity  of  the  heating  current  and  the  heating  time  (which 
changes  the  gauge  temperatures) .  Typical  data  from  these  tests  are  also 
plotted  in  Fig. 31  (labelled  as  wind-off  cases). 

It  was  found  that  the  heating  process  had  little  influence  on  the  temperature 
data  in  the  range  to  be  used  for  the  data  reduction  of  (3  -  7) .  This  view 
was  further  strengthened  by  measurements  of  the  heat  loss  (h^) .  For  all 
gauges  (whose  individual  resistances,  surface  areas  and  hence  heating 
conditions,  differed  widely)  it  was  found  that 

A 

—  h  -  0.1042  ±  0.0068  (3  -  8) 

s  w 

This  result  for  the  heat  loss  also  showed  that  the  thermal  capacity  of  the 
calorimeters  was  almost  in  proportion  to  the  gauge  surface  area,  since 
the  Pyrex  glass  plate  with  a  uniform  thickness  (0.058  ±  0.001  inches) 
was  expected  to  give  a  constant  heat  loss  over  the  test  surface.  This 
result  is  consistent  with  the  effective  thermal  capacity  discussed  in 
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Appendix  84.  The  equivalent  thermal  capacity  per  unit  area,  C  /A  , 

s  s 

2 

was  calculated  to  be  0.0414  Btu/ft  F  for  the  parabolic  thermal  model. 

Using  this  model  the  measured  heat  transfer  on  the  side-wall  without  the 
shock  (bare  plate  heat  transfer)  is  as  shown  in  Fig. 32.  In  the  figure 
the  local  Stanton  number  is  defined  as 

ch  -  V»i  “i  V  <3  -  9> 

The  average  local  Stanton  number  over  the  test  surface  is 

(C.  )  -  0.00188  ±  0.00027  (3  -10) 

hQ  avg. 

Although  the  measurements  show  a  15%  scatter  the  average  value  of 

-3 

1.88  x  10  agrees  well  with  the  prediction  from  the  Van  Driest  II 
theory  (Ref. 31)  of 

(C.  )  -  0.00181  (3  -11) 

hvD 

using  test  conditions  of 

Mach  number  ■  2.3  , 

stagnation  temperature  -  wall  temperature  »  285°K  , 

Reynolds  number  “  0.8  x  10^  and 
momentum  thickness  ■  0.040  inches. 

Using  this  flat-plate  boundary  layer,  the  influence  of  the  input  electric 
current  to  the  gauge  when  being  used  as  a  surface  thermometer  was  examined. 
These  data  are  also  shown  in  Fig. 31,  for  current  intensities  of  3  and  6mA. 
It  is  obvious  that  this  level  of  current  intensity  had  little 
influence  on  the  measured  data. 


-  35  - 


The  measurement  error  is  estimated  to  be  within  ±  15X.  This  figure  is 
a  little  worse  than  that  of  the  previous  slug-calorimeter  technique. 
However  it  must  be  remembered  that  the  present  stagnation  pressure  is 
only  about  £  of  that  used  in  the  intermittent  tunnel. 


3.6  Encapsulated  Liquid  Crystal  Technique 

To  visualize  the  global  temperature  distributions  on  the  side  wall  surface, 
encapsulated  liquid  crystals  were  used  in  the  manner  illustrated  in  Fig. 33. 
These  liquid  crystals  respond  to  small  differences  in  temperature  by 
showing  bright  colours.  When  illuminated  with  white  light  and  cooled 
through  their  sensitive  temperature  span,  they  show  a  continually  varying 
colour  spectrum  from  colourless  to  violet,  successively  through  blue, 
green,  yellow,  red  and  again  colourless.  This  colour  change  is  fully 
reversible.  The  model  painted  with  them  can  stand  repeated  use  without 
cleaning  and  recoating,  which  is  not  the  case  for  phase-change  paints. 

The  encapsulated  liquid  crystals  are  soluble  in  water  so  that  they  can  be 
sprayed  on  the  model  surface  to  form  a  layer  of  uniform  thickness.  Further 
information  can  be  found  in  Refs. 32  and  33* 

In  the  tests  conducted  in  the  intermittent  tunnel,  the  experimental 
procedure  was  very  simple.  Finished  matt-black  board  was  prepared  for  the 
side  wall.  The  encapsulated  liquid  crystals,  diluted  with  the  same  amount 
of  hot  water  (approximately  50°C),  were  sprayed  on  the  board  and  left  to 
dry.  Before  the  tunnel  run  the  board  was  always  allowed  to  reach  its 
equilibrium  temperature  (room  temperature).  When  the  tunnel  was  started 
and  the  air  began  to  cool  the  board  surface,  the  encapsulated  liquid  crystal 
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layer  showed  the  surface  temperature  patterns  changing  with  time.  The 
varying  coloured  pictures  of  the  surface  temperature  pattern  were  taken 
by  a  motor-driven  camera. 

The  sensitive  temperature  range  of  the  liquid  crystals  greatly  affected  the 
quality  of  the  pictures  obtained.  As  described  previously ,  this  tunnel 
provided  approximately  18°K  difference  between  the  wall  temperature  and  the 
adiabatic  wall  temperature.  The  sensitive  temperature  range  of  the  liquid 
crystals  used  was  chosen  so  as  to  have  a  5°K  span  at  approximately  5°K 
below  room  temperature;  i.e.  5-10°C  sensitive  liquid  crystals  for  the 
winter  season  and  18-23°C  sensitive  ones  for  the  sumner  season.  The  duration 
of  the  test  was  from  10  to  15  seconds  (for  the  hard-board  surface  used  above) . 

The  materials  of  the  side  wall  had  a  great  influence  on  the  quality  of  the 
pictures.  An  epoxy-resin  board  covered  with  copper  film  gave  uniform  thermal 
properties  over  the  test  surface.  However,  the  pictures  obtained  on  this 
board  tended  to  be  'smeared'  by  thermal  conduction  through  the  copper  film. 
Although  a  wooden  board  (birch)  gave  good  coloured  pictures,  fibres  of  the 
wood  surface  spoiled  the  surface  temperature  distributions.  Of  all  the 
surfaces  tested  it  was  found  that  the  smooth  face  of  standard  hard-board 
gave  the  best  results.  The  reference  surface  temperature  was  measured 
with  a  single  thin  film  gauge. 

The  liquid  crystal  layer  was  illuminated  at  approximately  45  degrees  to  avoid 
white-light  reflections  from  the  wall  surface  caused  by  the  small  grains 
of  the  encapsulated  liquid  crystals.  The  light  source  was  a  250  W  slide- 
projector.  The  camera  to  record  the  coloured  temperature  pattern  was  located 
so  as  to  take  pictures  nearly  normal  to  the  light  direction,  for  the  same 
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reason  as  above.  Kodak  Ektachrome  ASA  400  films  were  used  with  an  exposure 
time  of  1/60  seconds  at  an  aperture  of  f/5.6.  To  strengthen  the  colour 
contrast  of  the  pictures,  a  special  development  process  was  used  which 
increased  the  effective  film  speed  by  a  factor  of  4. 


3.7  Viscous  Layer  Surveys 

To  survey  the  boundary  layer  growing  along  the  tunnel  side  wall  and  the 
viscous  layers  in  the  interaction  region,  some  small  probes  and  a  probe 
drive  system  were  employed. 

The  probes  were  individually  designed  to  measure  the  total  pressure  (Pitot 
pressure),  the  static  pressure  and  the  flow  yaw-angle.  The  shapes  and 
dimensions  of  these  probes  are  shown  in  Fig. 34.  These  were  made  from 
hypodermic  tubing. 

The  probe  drive  system  was  designed  to  move  the  Pitot  probes  normal  to  the 
tunnel  side  wall  and  also  to  pitch  them  so  that  the  probes  could  be  aligned 
with  the  local  direction  of  the  flow.  Determination  of  the  distance  of  the 
probes  from  the  wall  surface  was  made  by  a  micrometer  built  into  the  drive 
system.  The  yaw  angle  of  the  probes  was  measured  directly  by  a  pointer 
moving  over  a  calibrated  scale. 


The  measured  total  and  static  pressures  were  used  to  obtain  other  properties 
in  the  viscous  layer.  The  method  of  analysis  is  described  in  Appendix  A. 
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4.  EXPERIMENTAL  RESULTS 

4.1  Presentation  of  Experimental  Results 

Some  typical  experimental  results  are  introduced  in  this  chapter,  together 
with  some  discussion  and  interpretation.  All  the  experimental  results  (for 
both  the  tunnel  test  programmes)  are  given  in  the  Appendices  as  listed  below. 

(1)  Schlieren  photographs  in  Appendix  C. 

(2)  Oil-flow  patterns  in  Appendix  D. 

(3)  Vapour-screen  pictures  in  Appendix  E. 

(4)  Wall-surface  static  pressure  distributions  in  Appendix  F. 

(5)  Wall-surface  heat  transfer  distributions  in  Appendix  G. 

(6)  Encapsulated  Liquid  Crystal  Pictures  in  Appendix  H. 

(7)  Wall-surface  temperature  distributions  in  Appendix  I. 

In  the  figures  two  different  types  of  Cartesian  co-ordinate  systems  are 

used,  of  which  the  origins  coincide  with  the  shock  generator  leading  edge 

on  the  tunnel  side  wall.  For  one  of  the  co-ordinate  systems,  the  X,  Y  and 

Z  axes  are  defined  in  directions  parallel  and  normal  to  the  free  stream. 

The  other  consists  of  the  X  and  Y  axes  which  are  respectively  defined 

8  8 

in  directions  along  and  normal  to  the  shock  generator  surface.  A  length 
Xg  is  also  employed  to  specify  a  distance  in  the  free  stream  direction 
measured  from  the  location  of  the  shock  wave.  The  details  are  shown  in  Fig. 35. 

4.2  The  Stability  and  Repeatability  of  the  Flow  Fields 

The  schlieren  system  was  continually  used  to  monitor  the  flow.  At  the 
beginning  of  each  test  programme  the  shock  generator  was  taken  through  the 
full  range  of  operating  angles.  Notes  were  made  of  any  unsteadiness 
resulting  from  either  shock  induced  separation  or  choking  of  the  flow. 
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For  shock  generator  angles  (4()  up  to  9°  with  Models  A  end  AM,  and  up  to 
15°  using  aodsl  B  the  flow  was  both  stable  and  repaatable.  A  comparison 
between  the  Measured  shock-wave  angles  and  those  predicted  by  the  inviscid, 
oblique- shock  theory  is  shown  in  rig  .36.  The  agreement  in  the  ranges  of  6^ 
given  above  is  good. 

With  model  A  the  flow  chokes  intermittently  at  6^  -  10°  and  chokes 
completely  when  6  •  11°.  For  model  AM  these  values  of  6g  are  lowered 
about  half  a  degree. 

For  model  B  inviscid  flow  theory  predicts  a  Mach  reflection  at  6g  -  16°  but 
this  was  not  observed  from  the  schlieren  picture.  However  the  oil-flow 
picture  at  <5g  -  15°  did  show  a  surface-flow  pattern  bending  round  to  form 
a  rather  unsteady  dead-flow  region  in  the  vicinity  of  the  bottom  wall.  This 
seemed  to  be  due  to  the  fully-deflected  side-wall  boundary  layer  interacting 
with  the  bottom  liner  (Appendix  D4-f ) . 

All  models  behaved  as  expected  in  the  useful  operating  range  of  6g  .  There 
were  no  spurious  compression  or  expansion  waves  to  upset  the  required  flow  field. 


4.3.  Flow  Field  Studies  Made  Using  the  Oil-flow  and  Vapour  Screen  Techniques 
4.3.1.  Oil-Flow  Patterns  on  the  Tunnel  Side  Wall 

As  the  shock  generator  angle  increased,  the  oil-flow  pictures  on  the  tunnel 
side  wall  showed  a  surface  flow  pattern  converging  toward  and  into  a  line 
near  the  shock.  Meanwhile  a  divergent  surface-flow  pattern  was  observed 
near  the  shock  generator  root,  for  all  shock  generator  angles.  For  these 
general  flow  features,  there  was  no  difference  between  Models  A,  AM  and  B. 


j 
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Typical  oil-flow  pictures  are  presented  in  Fig. 37,  together  with  sketches 
of  the  most  prominent  features. 

(1)  Surface-Flow  Convergence 

In  the  oil  flow  pictures,  two  types  of  surface-flow  convergence  are  observed. 
They  are  shown  in  Fig. 37,  by  using  solid  and  broken  lines  which  are  respectively 
named  complete  and  incomplete  convergence.  The  convergence  line  on  the 
side  wall  (whether  complete  or  incomplete)  is  labelled  (A)  in  the  sketches 
e.g.  Fig. 37. 

Complete  Surface-Flow  Convergence 

A  strong  surface-flow  convergence  into  a  single  line  appears  extensively  at 

6  -  9°  for  Models  A  and  AM  (e.g.  Figs. 37b  -  37c),  at  6  *  11°  for  Model  B 

®  8 

(e.g.  Fig.37f),  and  locally  near  the  shock  generator  at  6  <  9°  for  Models 

s 

A  and  AM  (e.g.  Fig. 37a) .  The  convergence  line  (indicated  by  a  solid  line) 
originates  in  the  vicinity  of  the  shock  generator  leading  edge  and  runs  a 
little  upstream  of  the  shock.  Ahead  of  this  convergence  line  the  oil  streak 
lines  completely  converge  into  it.  Downstream  of  the  line,  the  oil  coalescence 
is  not  clear  and  the  streak  lines  look  as  if  they  are  parallel  to  the  line. 
However,  in  this  region,  careful  observation  confirms  that  some  streak  lines 
do  in  fact  emerge  from  the  line  while  a  low  speed  surface-flow  area  is  formed 
between  the  convergence  line  and  the  shock  location  behind  it. 

Incomplete  Surface-Flow  Convergence 

In  the  picture  of  Model  B  at  6g  ■  11°  (Fig.37e),  it  can  be  seen  that  the 
•bove  complete  surface-flow  convergence  disappears  at  a  certain  distance 
from  the  shock  generator,  and  the  oil  streak  lines  near  the  convergence  line 
tend  to  swing  downstream  from  the  'disappearance  point'.  This  alteration 
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of  surface  flow  behaviour  along  the  convergence  line  can  be  found  in  most 
pictures  showing  some  complete  convergence,  (e.g.  Figs. 37b  -  37c  and  37e  -  37f). 
In  an  incomplete  convergence  region  the  oil  streak  lines  never  converge 
into  a  single  line,  but  they  still  tend  to  coalesce  near  the  shock  over  a 
'convergence  zone'  which  grows  in  size  along  the  shock.  The  centre  of  these 
streak  lines  converging  into  a  zone  is  indicated  by  a  broken  line  in  Fig. 37. 

The  surface-flow  pattern  forming  the  incomplete  convergence  zone  can  be 
seen  particularly  well  at  6g  *  10°  for  Model  B  (Fig.37d).  In  this  case,  the 
incomplete  convergence  line  originates  near  the  shock  generator  leading  edge. 

Thus,  there  is  a  limit  to  the  existence  of  complete  surface-flow  convergence 
in  terms  of  distance  from  the  shock  generator  leading  edge  and  this  distance 
depends  strongly  on  the  shock  generator  angle.  The  relation  between  this 
distance  and  the  shock  generator  angle  is  shown  in  Fig. 38. 

Boundary  Line  between  the  Inner  and  Outer  Surface  Flows 

The  convergence  line  described  above  (whether  complete  or  incomplete)  can 
be  used  as  a  boundary  line  to  divide  the  interaction  region  into  an  inner 
and  an  outer  surface-flow  region.  The  outer  surface  flow  is  upstream  of  the 
convergence  line  and  comes  from  the  upstream  undisturbed  surface  flow  of  the 
side-wall  boundary  layer. 

The  inner  surface  flow  originates  in  a  divergent  surface-flow  region 
(mentioned  below)  near  the  shock  generator  root  and  runs  downstream  of  the 
convergence  line.  It  is  clear  that  the  boundary  line  develops  from  the 
incomplete  convergence  line  to  the  complete  convergence  line,  as  the  shock 
generator  angle  is  increased. 
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(2)  Surface-Flow  Divergence 

In  the  oil-flow  pictures  on  the  side-wall,  a  surface-flow  divergence  is 
also  found  (at  any  shock  generator  angle)  along  the  shock  generator  root. 
The  divergent  surface-flow  pattern  is  typical  of  that  produced  by  'ordinary 
flow-reattachment'.  In  the  centre  of  the  divergent  surface-flow  region, 
it  is  possible  to  draw  an  'envelope'  of  the  oilstreak  lines.  In  Fig. 37 
the  line  is  marked  as  divergence  line.  The  divergence  line  starts  in  the 
vicinity  of  the  shock  generator  leading  edge  and  runs  downstream  whilst 
moving  away  slightly  from  the  shock  generator  root. 


(3)  Surface-Flow  Deflection  Angle 

As  the  shock  generator  angle  increases,  the  side-wall  surface-flow  deflection 
angle  grows  rapidly  and  eventually  exceeds  the  shock  angle,  as  observed  in 
the  oil-flow  pictures  already  discussed.  In  the  early  stage  of  the  growth 
(6g  >  9°  in  Fig. 37),  the  maximum  surface  deflection  is  close  to  that  of  the 
convergence  line.  As  the  surface-flow  convergence  line  strengthens 
to  form  a  complete  convergence  (e.g.  6g  >  9°  in  Fig. 37),  the  maximum  deflection 
angle  becomes  larger  than  the  convergence  line  angle.  In  this  case,  the 
maximum  flow  deflection  appears  near  the  shock  generator  root,  where  the 
above  surface-flow  divergence  is  observed.  The  surface-flow 

deflection  angle  is  plotted  in  Fig, 39,  with  respect  to  the  shock  generator 
angle,  together  with  the  surface-flow  angle  of  the  convergence  line.  In  this 
figure,  the  flow  angles  are  defined  relative  to  the  shock  generator  surface. 

The  figure  shows  that  the  maximum  surface-flow  deflection  angle  can  be  much 
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higher  than  McCabe's  prediction,  (see  Fig. 39).  It  must  also  be  stressed 
that,  even  when  the  surface-flow  deflection  exceeds  the  shock  angle  the 
surface-flow  pattern  does  not  always  show  the  complete  convergence  which 
is  one  of  the  features  needed  to  confirm  flow  separation.  In  other  words, 
the  condition  of  separation  or  even  'incipient  separation'  is  not  defined 
by  the  surface-flow  deflection  angle  exceeding  the  shock  angle. 


4.3.2.  Oil-Flow  Patterns  on  the  Shock  Generator  Surface 

In  Fig. 37,  the  oil-flow  pictures  on  the  shock  generator  surface  show  some 

interesting  features.  The  detailed  shock-generator  flow  patterns  for 

.  .  o 

Model  B  are  presented  in  Fig. 40.  Two  shock  generator  angles  (6  =10  and 

8 

13°)  are  shown  because  they  produce  the  incomplete  and  complete  types  of 
convergence  on  the  side-wall  surface  as  previously  described.  The  oil-flow 
pictures  (Figs.  37  and  40)  show  that  the  shock  generator  surface  flow  is 
divided,  by  two  narrow  compression  zones,  into  three  regions  of  convergent, 
accelerated  and  undisturbed  flow. 

The  oil-flow  pictures  (e.g.  Fig.  40a)  shows  that  a  narrow  compression  zone 
originates  from  the  shock  generator  leading  edge,  at  a  point  inside  the  side- 
wall  boundary-layer  edge,  and  spreads  downstream.  This  compression  zone 
divides  a  slender  flow-region  near  the  shock  generator  root.  In  this 
region,  a  convergent  surface-flow  is  observed,  and  is  labelled  convergence  B 
in  the  sketches.  This  convergence  region  grows  in  size  with  increasing 
shock  generator  angle.  When  the  side-wall  surface  flow  indicates  complete 
convergence  (convergence  A),  the  convergence  B  produces  a  narrow  foamy  zone 


developing  downstream.  This  convergence  B  is  caused  by  two  different 
surface  flows.  One  of  them,  which  is  inboard  of  the  convergence  line  B, 
is  connected  with  the  side-wall  surface  flow,  by  way  of  the  shock  generator 
root.  The  other  is  supplied  from  outboard  of  the  compression  zone. 

Outside  the  compression  zone,  an  accelerated  region  is  observed.  This  region 
lies  inboard  of  another  compression  zone  which  starts  from  the  leading 
edge  at  a  point  located  roughly  two  side-wall  boundary-layer  thicknesses 
away  from  the  shock  generator  root.  In  the  accelerated  region,  the  surface 
flow  is  accelerated  toward  the  shock  generator  root.  (The  outer  compression 
zone  is  not  clear  at  6^  -  13°  in  Fig. 40,  since  the  picture  was  taken  after 
the  tunnel  stop.  The  compressions  for  Models  A  and  AM  were  weak,  due  to  the 
existence  of  the  tunnel  liner  boundary  layer.) 

4.3.3.  Vapour-Screen  Pictures 

Typical  vapour-screen  pictures  are  presented  in  Fig. 41,  together  with 
sketches  showing  the  most  important  features.  In  the  figure,  the  approximate 
locations  of  the  shock  and  the  surface-flow  convergence  and  divergence 
regions  (as  found  from  the  corresponding  schlieren  and  oil-flow  pictures 
described  above)  are  all  marked. 

Shock  Wave 8 

In  the  vapour-screen  technique,  the  shock  is  rendered  visible  by  a  sudden 
change  of  the  fog  density  per  unit  volume  across  it,  the  number  of  fog 
particles  being  almost  proportional  to  the  air  density.  In  the  vapour-screen 
pictures,  the  shock  is  detected  as  a  boundary  line  between  lighter  and 
darker  regions,  as  shown  in  the  sketches.  The  shock  position  is  close  to 
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that  measured  from  the  corresponding  schlieren  picture,  in  the  inviscid 
external-flow  region.  The  shock  wave  approaches  the  side-wall  in  the  form 
of  a  lambda  shock  over  a  distance  from  two  to  four  side-wall  boundary-layer 
thicknesses. 

The  approximate  locations  of  the  shocks  as  measured  by  the  schlieren  and 
vapour-screen  techniques  are  compared  with  each  other,  in  Fig. 42.  The 
'vapour- screen  shock'  tends  to  appear  a  little  upstream  of  the  'optical  shock'; 
i.e.  the  vapour-screen  shock  angle  is  from  one  to  two  degrees  larger  than 
the  optical  one.  This  difference  corresponds  to  a  change  of  free-stream 
Mach  number  of  less  than  0.1  (for  ideal  air).  Although  the  reason  for  this 
difference  is  not  certain,  it  is  likely  that  the  moisture  added  to  the  tunnel 
air  for  the  vapour-screen  tests  changed  the  physical  properties  of  the 
test  flow. 

Small  Vortex  Swirling  near  the  Shock  Generator  Root 

The  vapour-screen  pictures  of  Model  B  (Figs. 41c  -  41e)  show  some  interesting 
features  near  the  shock  generator  root.  In  this  region  the  pictures  contain 
a  small  light,  circular  area.  Particularly  in  Fig. 41c,  the  light  circular 
area  can  be  seen  to  contain  a  darker  core.  (Although  the  other  pictures 
also  show  a  similar  region,  the  evidence  is  not  so  clear.)  The  light  circular 
area  looks  similar  to  those  produced  by  vortices  in  general  vapour-screen 
pictures,  (the  radial  acceleration  of  the  circular  flow  sweeps  away  the 
heavy  fog  particles  from  the  centre  of  the  flow  to  the  outer  region). 

On  the  side-wall  surface  'beneath'  the  small  light  circular  area,  there  is 
a  thin  dark  layer  spreading  across  the  wall.  This  shadowy  zone  is  similar 
to  regions  of  boundary  layer  flow  in  general  vapour-screen  pictures.  It  must 
be  emphasised  that  these  flow  features  are  observed  for  all  the  wedge 
generator  angles  tested. 
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With  the  help  of  the  corresponding  oil-flow  pictures,  the  following  inter¬ 
pretation  of  the  flow  field  is  made.  The  oil-flow  divergence  line  is 
located  near  the  junction  of  the  small  light  circular  area  with  the  thin 
dark  zone  (as  marked  in  the  vapour-screen  pictures) .  The  convergence  line  B 
is  near  the  other  side  of  the  light  circular  area,  on  the  shock  generator 
surface.  In  the  flow  region  near  the  shock  generator  root,  the  surface- 
flow  on  the  shock  generator  surface  (which  is  accelerated  toward  the  shock 
generator  root  and  converges  into  the  convergence  line  B)  lifts  up  from  the 
shock  generator  surface  at  the  convergence  line  B,  and  re-attaches  to  the 
side-wall  surface  at  the  divergence  line,  by  way  of  a  small  vortex  swirling 
near  the  shock  generator  root.  This  interpretation  also  explains  the 
appearance  of  divergent  flow  on  the  side-wall,  (which  exists  at  any  shock 
generator  angle,  even  when  no  surface-flow  convergence  is  observed  on  the  side 
wall) .  The  existence  of  a  small  vortex  near  the  shock  generator  root  would 
mean  that  the  divergence  line  is  the  re-attachment  line  'behind'  the  vortex. 

In  the  case  of  Model  AM  (Figs.41a-41b) ,  the  'thick  dark  zone'  comes  up 
onto  the  shock  generator  surface,  as  if  a  thick  'vortical'  flow  exists 
at  the  corner. 


Viscous  Layer 

As  noted,  the  viscous  layer  shows  as  a  thin  dark  zone  on  the  side-wall  surface 
It  is  further  observed  that  this  dark  zone  has  its  minimum  thickness  near 
the  location  of  the  convergence  line  A. 
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4.4.  Wall  Surface  Flow  Features  Indicated  by  the  Static  Pressure  and 
Heat-Transfer  Distributions 

4.4.1.  Static  Pressure  Fields 

The  isobar  patterns  on  the  side-wall  surface  are  presented  in  Fig. 43,  for 

Modal  A  at  6  ■  5,  7  and  9  degrees  and  for  Model  B  at  6  ■  7,  10  and 
s  s 

13  degrees.  These  isobars  were  obtained  by  interpolating  the  measured 
static  pressure  distributions  to  be  discussed  in  section  4.4.4. 

The  isobar  patterns  indicate  that  the  static  pressures  on  the  side-wall 
surface  gradually  increase  across  the  shock,  regardless  of  the  steep 
change  of  the  external-flow  pressure  expected  from  inviscid  theory. 

This  slow  growth  of  the  wall-surface  static  pressures  is  caused  by  the 
viscous  interaction  between  the  shock  and  the  boundary  layer  on  the  side- 
wall  surface. 

The  2.5  x  2.5  inch  wind  tunnel  test  programme  was  conducted  using  a 
relatively  small  test  flow  area.  From  the  results  a  suggestion  was  made  in 
Ref. 23  as  follows. 

The  static  pressure  on  the  side  wall  requires  approximately  ten  side-wall 
boundary- layer  thicknesses  to  reach  its  maximum  value,  and  the  upstream 
influence  due  to  pressure  propagation  through  the  subsonic  region  of  the 
viscous  layer  can  appear  up  to  approximately  six  side-wall  boundary-layer 
thicknesses  ahead  of  the  shock.  These  conclusions  were  made  on  the  basis 
of  the  quasi-two-dimensional  isobar  patterns  shown  for  example  in  Figs.  43a 
43c.  In  the  figures,  the  isobars  run  nearly  parallel  to  the  shock,  except 
in  a  region  near  the  tunnel  bottom  liner  where  the  shock  reflection 


disturbs  the  flow  field 
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The  next  test  programme  used  the  larger  test  flow  field  of  the  9x9  inch 
tunnel.  As  shown  in  Figs.  43d  -  43f,  the  isobars  in  the  larger  flow  field 
show  that  the  pressure  distribution  is  not  quasi-two-dimensional  but  more 
'conical'  in  nature.  In  a  small  limited  region  just  ahead  of  the  shock, 
the  isobars  run  roughly  parallel  as  in  a  quasi-two-dimensional  flow  but  overall 
it  can  clearly  be  seen  that  the  isobars  spread  over  the  side  wall  in  a 
roughly  conical  fashion  from  an  origin  near  the  shock  generator  leading  edge. 
Hence  it  is  not  particularly  useful  to  measure  the  pressure  propagation 
distance  in  terms  of  the  undisturbed  side-wall  boundary- layer  thickness, 
since  this  is  constant  whereas  the  pressure  field  is  growing  radially 
outwards  in  a  roughly  conical  fashion.  The  earlier  comments  concerning  the 
quasi-two-dimensional  nature  of  the  pressure  field  are  probably  conditioned 
to  the  small  test  flow  field  in  which  the  measurements  were  taken. 

(Similar  comments  can  be  made  for  the  oil-flow  patterns  in  Fig. 37. 

Although  the  complete  surface-flow  convergence  line  in  the  Model  A  conf iguration 
seems  to  run  parallel  to  the  shock,  it  behaves  in  a  more  'conical'  fashion 
for  Model  B  in  the  larger  test  flow  field.) 

It  must  be  stressed  that  the  basic  pattern  of  the  isobars  does  not  vary 
noticeably  with  increasing  shock  generator  angle.  For  both  Models  A  and  B, 
the  pattern  remains  the  same  but  the  pressure  gradient  across  the  shock 
is  steepened. 

4.4.2  Heat  Transfer  Distributions 

The  heat  transfer  patterns  over  the  test  surface  are  shown  in  Fig. 44,  where 
the  local  Stanton  number  is  presented  for  the  same  experimental  conditions 
as  for  the  pressure  fields  described  above.  These  patterns  were  again 
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obtained  by  interpolating  the  measured  heat  transfer  distributions  to  be 
discussed  in  Section  4.4.4. 

Fig. 44  gives  a  basic  pattern  of  the  heat  transfer  field  as  described  in  detail 
later.  The  basic  pattern  is  illustrated  in  Fig. 44a  on  which  the  most 
significant  surface-flow  features  from  the  oil-flow  pictures  are  superimposed. 
This  pattern  consists  of  three  individual  heat-transfer  peaks  which  can 
be  distinguished  from  each  other  as  follows. 

(i)  The  high  heat-transfer  'peak'  (A)  appears  near  or  a  little  upstream 
of  the  surface-flow  convergence  line. 

C ii)  Another  high  'peak'  ®  takes  place  closely  to  the  shock  generator  root, 
namely,  where  the  surface-flow  divergence  is  observed. 

(iii)  The  weak  'peak'  (5)  is  found  clearly  in  the  case  of  Model  A,  upstream 
of  the  peak  ® .  The  location  of  this  peak  is  a  little  upstream 
of  a  surface-flow  zone  where  the  oil  streak  lines  start  to  bend  to 
form  the  surface-flow  convergence. 

Model  A 

The  original  interpretation  (Ref .24)  of  the  heat  transfer  data  from  Model  A 

was  made  under  the  assumption  that  the  temperature  difference  (T  -  T  _) 

w  aw 

was  constant  over  the  interaction  region.  However,  more  recent  measurements 

for  Model  B  have  shown  that  it  is  necessary  to  take  account  of  the  local 

variation  of  (Tw  -  T#w),  particularly  downstream  of  the  shock,  as  discussed 

in  Section  5.1.  This  is  because  the  temperature  difference  was  small  for 

both  Models  A  and  B,  so  that  the  local  variation  of  T  (and  to  a  less  extent 

aw 

the  variation  in  T  )  cannot  be  ignored  in  analyzing  the  local  Stanton  number 
w 


from  measured  heat-transfer-rate  data.  It  is  now  believed  that  Figs. 44b  -  44d 
for  Model  A  underestimate  the  real  value  of  the  Stanton  number,  particularly 
downstream  of  the  shock. 

In  Figs.  44b  and  44c,  at  6  *  5  and  7  degrees  respectively,  the  heat 

8 

transfer  'peak'  (A)  appears  near  the  calculated  shock  location  near  where  the 
surface-flow  convergence  line  lies.  This  peak  is  represented  by  the  contour 
of  "  from  1.2  to  1.4  in  the  figures.  The  'peak'  (c)  is  found  upstream 

of  the  'peak'  (A)  as  separately  indicated  by  Cj/Sio*  1«2.  Downstream  of 
the  shock,  the  local  variation  of  (T  -  Tfiw)  gives  the  misleading  result 
of  the  local  Stanton  number  reducing  with  distances  from  the  shock.  This 
erroneous  behaviour  will  be  clarified  in  Section  5.1.  In  Fig.  44d  for 
6  *■  9°  ,  the  local  variation  of  (T  -  T  )  suggests  an  illusory  deep 

8  W  3W 

valley  between  the  shock  and  the  shock  generator  root.  However,  the  'peak'  (5) 
is  still  clearly  observed  close  to  the  shock  generator  root,  as  indicated 

*  VSo  *  1'4' 

Model  B 

The  local  Stanton  numbers  for  Model  B  were  analyzed  from  measured  heat-transfer 
coefficient  data.  (Refer  to  Appendix  B) .  In  this  method  of  analysis,  the 
local  variation  of  (T^  -  Tflw)  is  much  less  important.  Through  all  of  the 
shock  generator  angles  used,  ranging  from  5  -  7  to  13  degrees  (Figs.44e  -  44g) 

the  basic  pattern  of  the  heat  transfer  does  not  vary  noticeably  but  the 
heat  transfer  'peaks'  become  more  prominent  with  increasing  shock  generator 
angle. 

The  heat  transfer  'peak'  (X)  appears  just  ahead  of  the  surface-flow 
convergence  line,  as  represented  by  Ch/Cho>  1.0  in  Fig.44e  and  C^/C^  1.5 
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in  Figs.  44f  -  44g.  This  'peak'  (A)  is  located  more  upstream  of  the  shock 
and  covers  a  much  wider  area,  than  that  found  for  Model  A.  The  difference 
in  appearance  of  the  'peak'  (A)  between  Models  A  and  B  is  quite  consistent 
with  the  corresponding  oil-flow  patterns  (Fig. 37).  The  oil-flow  pictures 
show  that  the  surface-flow  convergence  in  the  test  conditions  of  Model  B 
takes  place  further  upstream  of  the  shock  and  is  less  concentrated  than  for 
Model  A.  It  is  concluded  that  the  appearance  of  'peak'  (X)is  closely  connected 
with  the  surface-flow  convergence. 

The  heat-transfer  'peak'  (S)  cannot  be  observed  in  Figs.44f  -  44g,  owing 
to  the  coarseness  of  the  'isotherms'  used.  However,  the  'peak*  (c)  can 
be  seen  from  the  heat  transfer  distributions  shown  in  Appendix  G,  Fig.G2. 

A  shallow  'trough'  appears  in  the  'peak  region'  (A)  as  if  it  separates 
the  'peak'  ©  from  the  main  'peak'  (K)  .  The  location  of  this  shallow 
'trough'  approximately  coincides  with  the  surface  stream  lines  beginning 
to  bend  towards  the  plane  of  the  shock.  This  heat  transfer  behaviour  is 
almost  the  same  as  in  Model  A  showing  clear  distinguishment  between  the 
'peaks'  (A)  and  (c)  . 

The  high  heat  transfer  'peak'  (B)  is  very  evident  in  Figs.44f  -  44g  with 
Cl/Cho  reach^n8  values  as  high  as  3.  The  heat  transfer  in  this  region  close 
to  the  shock  generator  root  will  be  discussed  in  detail  in  Section  4.4.4. 


4.4.3.  Surface  Temperatures 

The  surface  temperature  patterns  are  presented  in  two  ways;  i.e.  by  using 
the  encapsulated  liquid  crystal  pictures  for  Model  A,  and  surface  temperatures 
measured  with  the  thin-film  gauges  for  Model  B. 
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Typical  encapsulated  liquid  crystal  pictures  are  shown  in  Fig. 45,  for 
6g  ■  7  and  9  degrees.  In  the  figure,  the  'blue-tone'  and  'yellow-tone'  zones 
shows  higher  and  lower  temperature  regions  respectively.  (In  Fig. 45,  the 
'blue-tone'  and  'yellow-tone'  zones  are  respectively  shown  as  'dark'  and 
'light'  greys.  Since  the  colour  photographs  show  only  two  main  colours 
they  have  been  printed  in  black  and  white  for  convenience.)  The  colour 
patterns  change  with  time  after  the  tunnel  starts.  The  time  indicated  in 
the  figure  is  within  the  period  when  the  colour  pictures  give  the  most 
information  about  the  general  surface-temperature  patterns.  The  detailed 
temperature-time  histories  are  presented  in  Appendix  H. 

It  is  suggested  that  the  higher  and  lower  temperature  zones  correspond  to 
the  lower  and  higher  heat-transfer-rate  regions  respectively.  This  is 
because  the  temperature  patterns  were  obtained  under  a  'hot  wall'  condition 
at  the  moment  that  the  intermittent  tunnel  started,  and  the  thermal  capacity 
of  the  wall-surface  material  is  regarded  as  constant  over  the  test  surface. 
Then,  if  the  thermal  conduction  on  the  wall  surface  can  be  ignored,  the 
wall  surface  is  cooled  in  proportion  to  the  product  of  heat  transfer 
coefficient  (h)  and  temperature  difference  (Tw  -  Tgw> .  This  idea  makes  it 
possible  to  compare  these  surface  temperature  patterns  with  the  heat  transfer 
patterns  based  on  heat- transfer-rate  data  (Figs. 44b  -  44d) . 

One  can  clearly  find  some  correspondence  between  the  colour  and  heat  transfer 
patterns.  The  colour  patterns  show  lower  temperature  regions  at  two  locations 
for  all  the  shock  generator  angles  tested;  i.e.  near  the  shock  and  near  the 
shock  generator  root,  (Fig. 45).  One  of  the  lower  temperature  regions 
indicates  that  a  high  heat  transfer  zone  exists  near  the  shock.  This  high 
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heat-transfer  zone  corresponds  to  the  heat  transfer  'peak'  ®  (Fig. 44a). 

The  encapsulated  liquid  crystal  pictures  also  suggest  that  the  'peak'  (A) 
could  be  spreading  more  extensively  along  the  shock  than  in  Figs. 44b  -  44d. 

The  second  lower  temperature  zone  corresponds  to  the  high  heat-transfer 
'peak'  (§)  lying  close  to  the  shock  generator  root.  However,  it  must  be 
stressed  that  these  coloured  surface-temperature  patterns  cannot  directly 
give  absolute  quantity  of  the  local  Stanton  number  since  they  are  also 
influenced  by  the  local  variations  of  the  temperature  difference  (Tw  -  T  w) . 

The  measured  surface-temperature  patterns  for  Model  B  are  presented  in  Fig. 46. 
These  patterns  indicate  a  low  temperature  zone  between  the  shock  and  the 
surface-flow  convergence  line,  for  all  shock  generator  angles,  as  indicated 
by  Twl  -  T^  <  -  2°K.  A  high  temperature  region  is  observed  downstream 
of  the  shock,  as  represented  by  Twl~  T^  >  +  2°K.  These  temperature  data 
will  be  used  in  Section  5.1.  to  discuss  the  heat  transfer  data  measured 
across  the  shock. 


4.4.4.  Static-Pressure  and  Heat-Transfer  Distributions 

Typical  wall-surface  static  pressure  distributions  in  the  free-stream  or 

external-flow  directions  are  presented  in  Fig. 47. 

The  pressure  distributions  plotted  in  Fig. 47a  for  Model  A,  were  measured 
in  the  free  stream  direction  on  the  tunnel  centre  line  (Y  -  1.25  inches). 
These  distributions  show  that  the  wall-surface  static  pressure  gradually 
increases  across  the  shock.  No  noticeable  alteration  of  the  basic 
distribution  pattern  is  observed  for  shock  generator  angles  ranging  from 
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5  to  9  degrees,  though  the  corresponding  oil-flow  pictures  indicate 
considerable  change  of  the  surface-flow  pattern  from  the  incomplete 
convergence  to  the  complete  convergence  (Figs. 37a  -  37b). 

The  pressure  distributions  in  Fig. 47b  for  Model  B  were  measured  along  lines 

parallel  to  the  shock  generator  surface.  In  this  figure,  a  slight  alteration 

of  the  distribution  pattern  is  observed  as  the  shock  generator  angle  is 

changed  from  10  to  13  degrees.  The  region  near  the  shock  position  tends 

to  be  'flattened'  with  increasing  shock  generator  angle.  This  tendency 

appears  more  strongly  at  Y  >2.04  than  Y  >1.04  inches.  The  alteration 

8  8 

of  the  pressure  distribution  pattern  is  consistent  with  the  change  of  oil-flow 
pattern  (Figs.37d  -  37f ) .  As  the  shock  generator  angle  increases  from 
10  to  13  degrees,  the  incomplete  convergence  of  the  surface-flow  changes 
to  a  complete  convergence,  and,  with  increasing  distance  from  the  shock 
generator,  the  complete  convergence  line  is  followed  by  an  extensive  low-speed 
zone  between  itself  and  the  shock  located  downstream  of  it.  Fig. 47b  shows 
that  this  is  precisely  the  region  in  which  the  pressure  is  'flattened' 

(rising  slowly).  The  positions  of  the  complete  convergence  line  and  the 
shock  are  marked  on  the  figures.  The  surface  pressure  only  begins  to  rise 
rapidly  again  once  it  has  passed  the  shock  location. 

In  Fig. 48,  typical  heat-transfer  and  static-pressure  distributions  in  the 
direction  normal  to  the  free-stream  or  the  external-flow  are  presented. 

Locations  of  the  shock  and  the  surface-flow  convergence  line  are  indicated  on 
the  figure.  As  described  previously,  the  static  pressures  gradually  increase 
across  the  shock.  However,  the  heat  transfer  distributions  have  very  complicated 
patterns.  (In  Figs. 48a  -  48b  for  Model  A,  a  high  heat  transfer  region  is 
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This  flattening  of  the  pressure  distribution  is  thought  to  be 
an  indication  of  sudden  viscous  layer  thickening. 
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observed  near  the  tunnel  bottom  liner  at  Y  -  2.5  inches.  Since  this  is 
caused  by  interference  from  the  tunnel  liner,  the  present  discussion  omits 
this  region.) 

Surface-Flow  Convergence  Region 

The  surface-flow  convergence  brings  about  the  high  heat  transfer  'peak'  (X) 
near  itself  or,  in  most  of  cases,  outside  it.  However,  no  noticeable  peak 
of  the  static  pressures  is  detected  in  this  region.  Although  Scuderi  (Ref .34) 
proposed  a  power  law  correlation  between  the  heat  transfer  and  the  static 
pressure  in  this  region,  the  present  experimental  data  do  not  show  any 
correlation  between  them  at  all.  As  noted  in  Section  4.4.2.,  the  heat- 
transfer  'peak'  (A)  is  connected  only  with  the  surface-flow  patterns. 

Interaction  Region  downstream  of  the  Shock 

Downstream  of  the  shock,  >oth  the  static  pressure  and  the  heat  transfer 
tend  to  increase  as  the  shock  generator  surface  is  approached,  (except  in 
Fig. 48b  where  the  heat  transfer  is  considered  to  be  under-estimated  because 
of  the  local  variation  of  (T  -  T  ),  as  noted  earlier).  Neumann  et  al. 

(e.g.  Ref. 17)  proposed  a  correlation  between  the  static  pressure  and  the 
heat  transfer  in  this  region.  This  correlation  is  tried  in  Fig. 49. 

In  Fig. 49a,  the  local  Stanton  numbers  are  plotted  with  respect  to  the  static 
pressure  and  compared  with  the  well-known  power  law  correlation  for  two- 
dimensional  aerodynamic  peak  heating: 
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h  max 
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where  n  ■  0.8  for  a  turbulent  boundary  layer.  (Refer  to  Back  et  al. ,  Ref. 35 
or  Neumann,  Ref.l)  .  However,  the  present  experimental  results  show  little 
correlation  between  pressure  and  heat  transfer  in  this  region.  In  Fig. 49b, 
the  maximum  heat  transfer  of  the  'peak'  (b)  is  considered  in  streamwise 

direction,  while  correlated  with  the  corresponding  static  pressure. 

/0.8 

'P/P,)  *  }  and  n  are  plotted  versus  non-dimensional 
l  max 

streamwise  distance  using  the  undisturbed  side-wall-boundary-layer  thickness. 
The  figure  shows  that  the  three-dimensional  interaction  can  give  rise  to 
higher  heat  transfer  than  the  two-dimensional  case,  particularly  near  the 
shock  generator  leading  edge  region. 

The  heat  transfer  in  the  three-dimensional  interaction  region  will  be  further 
discussed  in  Section  5. 


4.5  Viscous  Layer  Surveys 

Typical  viscous-layer  properties  measured  using  static  and  pitot  pressure 
probes,  and  a  yaw  meter,  are  presented  in  Fig. 50  for  Model  A  with 
&a  *  5,  7  and  9  degrees.  The  experimental  results  are  presented  using  the 
(Xg»  Yg.  Z)  co-ordinate  system.  Thus  the  streamwise  and  cross-flow  velocity 
components,  u  and  v  respectively,  are  defined  in  directions  parallel 
and  normal  to  the  external  flow  downstream  of  the  shock.  The  survey  results 
show  that  the  side-wall  boundary  layer  is  skewed  a  great  deal  and  warped  near 
the  shock.  The  salient  features  of  the  flow  are  as  follows. 
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(1)  The  undisturbed  side-wall  boundary  layer  upstream  of  the  shock,  has  no 
detectable  flow-deflection  and  shows  a  velocity  profile  typical  of  & 
fully  turbulent  boundary  layer.  Thus,  as  expected,  the  flow  is  quasi- 
two-dimensional.  (Refer  to  Fig. 16).  In  the  interaction  region  down¬ 
stream  of  the  shock,  this  boundary  layer  has  been  changed  dramatically 
into  a  complex  three-dimensional  zone  of  highly  skewed  flow.  The  flow 
deflection  angle  is  zero  in  the  inviscid  external  region,  because  the 

external  flow  is  parallel  to  the  shock  generator  direction  (x  axis) . 

8 

As  the  side  wall  is  approached  so  the  flow  direction  yaws  or  skews. 

(2)  At  6^  ■  5  and  7  degrees,  the  local  flow  deflection  in  the  outer  part  of 

the  boundary  layer  inclines  lightly  towards  the  shock  generator,  as  the 

distance  from  the  wall  is  decreased.  (See  0.2  inches  <  Z  <  0.4  inches 

in  Fig. 50a).  In  this  flow  region,  the  local  static  pressure  changes 

down  rapidly  from  the  external-flow  pressure  to  the  wall-surface  pressure 

(Fig. 50b)  and  the  local  Mach  number  exceeds  the  external-flow  Mach 

number  (Fig. 50c).  This  flow  behaviour  suggests  that  the  three-dimensional 

glancing  shock  interaction  provides  a  flow  region  similar  to  the 

'supersonic  tongue'  of  two-dimensional  normal  shock  interaction. 

(Refer  to  Fig. 7).  In  the  case  6  ■  9°  a  similar  flow  behaviour  can  be 

8 

detected  (weakly)  in  the  region  of  0.3  inches  <2  <  0.6  inches. 

(3)  In  the  inner  region,  2  <  0.2  inches  for  »  5  and  7  degrees,  the  local 
flow  deflection  angle  increases  steadily  as  the  wall  is  approached  and 
reaches  the  surface-flow  deflection  angle  (obtained  from  the  oil-flow 
picture)  at  2  -  o,  (Fig. 50a).  In  this  flow  region,  the  local  static 
pressure  is  eloost  constant  (Fig. 50b)  while  the  local  Mach  number  simply 


decreases  to  zero  at  the  wall  (Fig. 50c).  The  corresponding  velocity 
profiles  show  the  typical  patterns  of  a  highly  skewed  'boundary  layer* 
(Figs. 50d  -  50e) .  The  cross-flow  velocity  component  is  considerable 
near  the  wall. 

(4)  The  distribution  of  the  local  flow  deflection  angle  for  6s  ■  9°  is  in 

striking  contrast  to  those  at  6  ■  5  or  7  degrees  (Fig. 50a). 

8 

Below  Z  ■  0.1  inches,  the  local  deflection  angle  increases  rapidly, 
overshoots  the  oil  flow  angle  at  the  surface,  then  swings  back  to 
match  it  at  Z  •  0.  The  corresponding  streamwise  velocity  profile  has 
a  point  of  inflection  and  a  lower  region  of  reduced  velocity  (Fig.50d). 


I 
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5.  HEAT  TRANSFER  IN  THE  INTERACTION  REGION 

5.1.  Heat  Transfer  Measurement  in  the  Shock  Region 

In  the  previous  section,  it  was  mentioned  that  local  variations  of  the 
temperature  differences  (T  -  Taw) »  i.e.  between  the  wall  and  adiabatic 
wall  temperatures,  greatly  affected  the  interpretation  of  measured  heat- 
transfer-rate  data,  particularly  when  the  absolute  value  of  (T  -  T&w)  is 
small.  This  is  because  the  local  Stanton  number,  C^,  is  defined  by 

ch  ■  i/piuicP  <T„  -  T.„>  «  -  5> 


where  q  is  the  measured  heat  transfer  rate  while  p,,  u,  and  C  are  the 

11  p 

free  stream  density,  the  free  stream  velocity  and  the  specific  heat 
respectively.  However,  the  local  variation  of  (Tw  -  Tflw)  is  much  less 
important  if  the  local  Stanton  number  is  obtained  from  measured  heat-transfer- 
coefficient  data,  h  ,  as  given  by 


C 

h 


h./piuicP 


(3  -  9) 


which  does  not  involve  the  term  (T  -  T  ) .  C.  was  derived  from  (3  -  5)  in 

w  aw  h 

the  Model  A  tests  but  from  <3  -  9)  in  the  Model  B  tests. 


During  the  heat  transfer  measurements  on  Model  B,  the  local  wall-surface 
temperature ,  T^,  was  also  measured  (Fig. 46).  This  temperature  data 

gave  the  local  variation  of  (Tw  -  Taw) ,  i.e.  (Twl  -  T  j) ,  through  a 
thermal  balance  on  the  wall  surface  (Appendix  B4) ;  i.e. 


h.(I«l  -  W  ■  VW 


(the  heat  transfer  from  (The  heat  transfer 
the  calorimeter  to  the  from  the  inner  wall  to 
air)  the  calorimeter) 


(5  -  1) 
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where  is  the  reference  wall  temperature  and  h_  is  a  heat-loss  coefficient 
The  results  of  (T^  -  Tawj)  are  presented  in  Fig. 51,  for  the  bare-plate  case 
and  the  interaction  cases  at  6  =  7,  10  and  13  degrees.  (In  the  present 
results,  the  measured  wall  reference  temperature  is  adjusted  by  using  the 
'undisturbed'  free-stream  adiabatic-wall  temperature  assuming  a  recovery 
factor  of  r  *  0.89).  In  the  figure  (T^  -  T  p  for  the  bare-plate  case 
is  given  by  the  solid-round  symbols  and  shows  an  almost  flat  distribution 
over  the  test  surface.  However,  for  the  interaction  cases,  (T  -  T  ,) 
has  a  significant  local  variation  downstream  of  the  shock;  i.e.  more  than 
10°K.  This  local  variation  is  naich  larger  than  the  expected  change  of 


local  adiabatic-wall  temperature  across  the  shock,  which  is  estimated 

to  be  2  -  4°K  for  6  *  7°  -  13°  assuming  a  constant  recovery  factor  of 

s 


r  *  0.89.  The  actual  recovery  factor  has  been  calculated  from  the  measured 
temperature  data  by 


T  -  1/(1  +  M2) 

«  o  z 


(5  -  2a) 

(5  -  2b) 


where  Tq  is  the  stagnation  temperature  and  y  the  ratio  of  specific  heats, 
while  M  is  taken  to  be  local  external-flow  Mach  number;  i.e.  and 
upstream  and  downstream  of  the  shock  respectively.  The  recovery  factor 
distributions  are  presented  in  Fig. 52a,  and  the  patterns  over  the  test 
surface  are  shown  in  Figs. 52b  -  52d.  Although  the  accumulated  calculation 
errors  reduce  the  accuracy  of  the  results,  the  recovery-factor  behaviour  is 
fairly  clear.  As  the  flow  goes  downstream  the  recovery  factor  decreases 
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in  the  region  between  the  surface-flow  convergence  line  and  the  shock,  but, 
downstream  of  the  shock,  it  quickly  increases  almost  in  proportion  to  the 
distance  from  the  shock  and  approaches  unity. 

It  can  now  be  estimated  how  much  the  local  variation  of  (T  -  T  _)  affects 

w  aw 

the  local  Stanton  number  previously  obtained  from  the  heat-transfer-rate 
data  of  Model  A  by  assuming  that  (T^  -  is  equal  to  constant. 

From  (3  -  5)  and  (3  -  9)  it  can  be  shown  bhat 


°h 


(5  -  3) 


where  is  the  corrected  local  Stanton  number  and  is  the  value  computed 
assuming  (Tw  -  Tfiw)  to  be  constant.  When  the  local  value  of  the  recovery 
factor  is  known,  the  variation  of  (Tw  -  Tftw)  can  be  calculated  from 


<T 


wl 


-  T  ) 
awl  cal. 


(T  -  T  ) 
w  o 


+  (1  -  r)(To-  T J 


(5  -  4) 


If  the  absolute  value  of  (T^  -  Tq)  is  large  enough  to  neglect  the  second 
term  in  (5-4),  then  it  is  clear  from  (5  -  3)  that  can  be  regarded 
as  equal  to  C^.  However,  in  the  present  test  conditions  for  Model  A, 

(T  -  T  )  *  5°K  at  the  maximum,  so  the  second  term  cannot  be  ignored. 

w  O 

This  second  term  can  be  estimated  to  vary  locally  from  7.5°K  to  22.5°K 
for  the  recovery  factor  ranging  from  0.85  to  0.95  (as  expected  from  Fig. 52 
for  Model  B) ,  and  give  approximately  ±  40Z  ill-interpretation  of  the  local 
Stanton  number  under  the  assumption  of  constant  (Tw  -  . 
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5.2.  Typical  Heat  Transfer  Distributions 

Some  typical  heat  transfer  results  for  Model  B  are  plotted  in  Fig. 53  and 
compared  with  Neumann's  results  (Ref. 17),  which  were  measured  at  almost 
the  same  shock  pressure  ratio. 

In  the  interaction  region  between  the  shock  and  the  shock  generator  root, 
there  is  good  agreement  between  the  three  sets  of  data.  In  all  of  the  cases, 
the  heat  transfer  rises  to  'peak'  (§) near  the  shock  generator  root.  However, 
in  the  region  ahead  of  the  shock  (y/Y  >1  in  Fig. 53),  the  agreement  is 
confused  by  various  patterns  of  the  'peak'®.  The  heat-transfer  results 
in  this  region  probably  indicates  how  important  the  individual  test  conditions 
are  and  perhaps  how  only  the  grossest  flow  characteristics  are  likely  to 
be  coimnon  when  comparing  different  test  data  in  different  wind  tunnels. 
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THREE-DIMENSIONAL  FLOW  SEPARATION  IN  THE  INTERACTION  REGION 

6.1.  Definition  of  Three-Dimensional  Separation 

As  noted  previously,  the  oil-flow  pictures  on  the  side  wall  show  that  the 
surface  flow  in  the  interaction  region  is  highly  deflected  with  increasing 
shock  generator  angle  and  forms  a  convergence  zone  near  the  shock  (convergence 
A  in  Fig. 37).  When  the  shock  generator  angle  exceeds  a  certain  value,  this 
surface-flow  convergence  zone  produces  a  'complete  convergence  line'  a 
little  upstream  of  the  shock.  See  Fig. 54.  In  front  of  this  line,  the  oil- 
streak  lines  coalesce  with  it.  However,  behind  the  line,  the  oil-streak 
lines  run  nearly  parallel  to  it.  This  oil-flow  behaviour  is  as  shown  by 
'convergence  pattern  I'  in  the  figure.  This  is  not  like  'convergence 
pattern  II'  indicating  the  surface-flow  coalescence  into  a  single  line  on 
both  sides  of  it,  which  is  typical  of  Maskell's  ordinary  flow  separation 
(Fig. 6b,  taken  from  Ref. 6)  and  commonly  accepted  to  be  a  type  of  three- 
dimensional  separation.  However,  the  present  paper  shows  that  the  surface- 
flow  convergence  pattern  I  also  belongs  to  the  same  category  of  ordinary 
flow  separations  as  Maskell's.  This  assertion  is  based  on  Lighthill’s 
'skin-friction-line  topology'  (see  Ref .37  or  the  excellent  review  of 
Tobak  et  al.,  Ref. 38),  as  follows. 

The  flow  separation  is  easily  and  clearly  defined  by  using  the  'limiting- 
stream-line  concept'.  The  limiting  stream  line  is  the  closest  stream  line 
to  the  wall  surface  in  the  attached  flow  field.  Thus  for  attached  flow  this 
stream  line  coincides  with  the  'surface  stream  line’  (which  is  employed  in 
other  parts  of  this  paper)  and  is  in  contact  with  the  'skin-friction  line' 
(surface  shear-stress  trajectory).  Three  dimensional  separation  can  be 
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defined  as  follows.  When  the  limiting  stream  line  leaves  the  wall  surface, 
the  main  flow  separates  from  the  wall  at  the  separation  point  of  the 
limiting  stream  line.  Now  having  given  a  clear  definition,  what  is  required 
is  proof. 

The  limiting-stream-line  concept  is  not  sufficient  to  give  conclusive  evidence 
of  three-dimensional  flow  separation  in  the  oil-flow  pictures,  for  the 
following  reasons. 

After  the  limiting  stream  line  rises  up  from  the  wall  surface,  no  trace  of 
its  behaviour  can  be  obtained  on  the  wall  surface.  The  separation  line 
lies  exactly  on  a  line  at  which  the  limiting  stream  line  disappears  from  the 
wall  surface.  Even  after  the  disappearance  of  the  limiting  stream  line 
from  the  wall  surface,  the  oil-streak  line  is  left  on  the  surface  and  moves 
with  the  skin  friction.  Thus,  the  oil-streak  lines  cannot  directly  indicate 
behaviour  of  the  limiting  stream  line  in  the  vicinity  of  the  separation  line. 

Introduction  of  the  ' skin-friction-line '  concept  makes  it  possible  to 
clarify  the  flow  behaviour  in  the  vicinity  of  the  separation  line.  The  skin- 
friction  line  is  defined  everywhere  on  the  wall  surface  and  the  separation 
line  coincides  with  a  skin-friction  line  (showing  an  envelope  of  an  infinite 
number  of  skin-friction  lines).  Using  Lighthill's  simple  argument  it  is 
possible  to  discuss  how  the  limiting  stream  line  leaves  the  wall  surface  in 
the  vicinity  of  the  separation  line.  A  stream  tube  is  supposed  to  be  on 
the  wall  surface  near  the  separation  line  as  shown  in  Fig. 55a.  The  width  and 
height  of  the  rectangular  stream  tube  are  n  and  h  respectively  (which  are 
assumed  to  be  finite  but  nearly  'infinitesimal'  so  that  the  local  resultant 
velocity  vectors  are  coplanar  and  form  a  linear  profile) .  Then  the  mass 
flux  through  the  stream  tube  is  given  by 

i 
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m  “  phnU 


(6  -  1) 


where  p  is  the  density  and  U  the  mean  resultant  velocity  in  the  cross 
section  of  the  stream  tube.  The  skin-friction  may  be  written  as 

* 


T 


,U  \ 


(6  -  2) 


where  p  is  the  coefficient  of  viscosity  so  that 
U*  -  th/2p 


(6  -  3) 


Hence, 


m 


ph  nt/2 p 


constant 


yielding 


(6  -  4) 


(6  -  5) 


Equation  (6  -  5)  shows  that  the  height  of  the  limiting  stream  tube  on  the 
wall  increases  rapidly  towards  infinity,  when  either  the  width  (n)  or  the 
skin-friction  (t)  approaches  zero.  Therefore,  if  any  limiting  stream  line 
is  situated  in  the  stream  tube,  the  limiting  stream  line  must  leave  the  wall 
surface  (and  the  flow  is  separated)  in  either  of  the  following  cases: 


(i)  The  skin-friction  vanishes  on  the  wall,  as  in  Maskell's  singular 
separation  at  a  singular  point,  (Fig. 6a). 

(ii)  The  width  of  the  stream  tube,  namely  the  distance  between  the 
adjacent  friction  lines  becomes  zero,  as  in  Maskell's  ordinary 
separation  along  the  separation  line  (Fig. 6b). 


The  latter  case  clearly  occurs  in  the  surface-flow  convergence  pattern  I  in 
Fig. 54,  upstream  of  the  complete  convergence  line.  The  distance  between  adjacent 
oil-streak  lines  indicating  the  skin-friction  lines  reaches  zero. 
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The  side-wall  boundary  layer  upstream  of  the  complete  convergence  line 
must  be  separated  at  the  line,  without  the  need  to  refer  to  the  downstream 
surface-flow  pattern. 


6.2.  Angular  Momentum  Conservation  Concept  of  the  Surface  Flow  Convergence 
6.2.1.  Consideration  of  the  Oil-Flow  Convergence 

The  oil-flow  pictures  on  the  side  wall  show  the  following  behaviour  of 
the  surface-flow  convergence  near  the  shock  (convergence  A) ,  as  described 
in  Section  A. 3.1. 

(i)  The  surface  flow  can  preserve  an  incomplete  convergence  pattern 
(showing  attached  flow),  even  if  the  surface-flow  deflection  exceeds 
the  shock  angle  (Fig.37d). 

(ii)  Even  when  the  shock  generator  angle  is  large  enough  to  produce  surface- 
flow  convergence  (showing  separated  flow),  the  complete  convergence 
line  can  change  at  a  certain  distance  from  the  shock  generator  into 

an  incomplete  convergence  zone  which  bends  round  in  the  downstream 
direction.  (Figs. 37b  -  37c  and  37e  -  37f). 

These  experimental  facts  suggest  that  the  appearance  of  flow  separation 

in  the  interaction  region  does  not  depend  on  the  surface-flow  deflection 

angle  exceeding  the  shock  angle  (as  proposed  by  McCabe,  Ref. 8),  but  on  other 

flow  conditions.  The  suggestion  made  here  is  that  there  is  a  limit  to  the 

angle  through  which  the  lower  part  of  the  boundary  layer  can  be  deflected 

or  skewed.  The  maximum  permissible  angle  of  the  surface-flow  deflection 

is  being  discussed  in  the  following  Sections. 


6.2.2.  Review  of  McCabe's  Theory 

McCabe  (Ref. 8)  developed  a  simple  approximate  theory  which  related  the 
surface-flow  deflection  to  the  external  flow  conditions.  In  this  theory, 
circulation  around  a  rectangular  vortex  tube  is  considered  in  a  quasi-two- 
dimensional  flow  field  assuming  that  the  boundary  layer  thickness  is 
infinitesimal  so  that  the  flow  can  slip  at  the  wall  surface  (Fig. 56a).  The 
vortex  tube  has  a  height  equal  to  the  boundary- layer  thickness  and  is 
initially  situated  upstream  of  the  shock  so  as  to  be  normal  to  the  free 
stream  in  a  spanwise  sense.  Then,  the  circulation  around  this  vortex  tube 
is  given  by 

(li  6xidz  "  ui6xi  (6 " 6) 

boundary  boundary 

layer  layer 

where  da^  is  local  vorticity  for  a  vortex-tube  element  with  height  dz^ 
(normal  to  the  wall)  and  width  dx^  (parallel  to  the  free  stream) ,  while  z 
is  the  displacement  from  the  wall,  u  the  streamwise  velocity  and  u^  the 
free-stream  velocity.  This  vortex  tube  travels  downstream  at  the  external 
flow  velocity  and  passes  through  the  shock  into  the  interaction  region. 

In  this  process,  the  vortex  tube  is  assumed  to  retain  a  rectangular  cross 
section  with  its  height  constant  but  its  width  changing  to  maintain 
continuity  of  the  flow.  Downstream  of  the  shock,  the  vortex  tube  axis  has 
been  swept  relative  to  the  external  flow  direction,  owing  to  the  existence 
of  oblique  shock,  so  unit  spanwise-length  of  the  vortex  tube  has  been  also 
changed,  (Fig. 56a).  In  this  flow  region,  the  circulation  around  the  vortex 
tube  is  given  in  the  external  flow  direction,  as 


-  68  - 


I  Sw.cosc  6x«tanC  6  *  l  -  t—  6x_tan£  dx  ■  v  dx.tanC  (6  -  7) 

L  l  Z  )  oZ  l  W  L 

boundary  boundary 

layer  layer 

where  du^,  dxj  and  Sz2  are  downstream  value  corresponding  to  dw^,  dx^ 
and  dz^  mentioned  above,  c  is  the  sweep  angle  o£  the  vortex  tube  relative 
to  the  external  flow  and  is  the  cross-flow  velocity  at  the  wall. 
Conservation  of  the  circulation  needs  (6  -  7)  to  be  equal  to  (6  -  6)  ,  then 

Idw^dx^dz^  ■  1 60)26x28105  6x2  (6  “  8) 

boundary  boundary 

layer  layer 

so  that 

UjdXj^  ■  v^dxgtanc  (6  -  9). 


The  surf ace-flow  deflection  angle  relative  to  the  external  flow,  c,  is 
defined  by 

tane  ■  v  /u„  (6  -10) 

w  2 

where  U2  is  the  downstream  external-flow  velocity. 

Substituting  (6  -  9)  into  (6  -10)  and  calculating  the  changes  of  flow 
properties  across  the  shock,  gives 

cos^(8-6  )-  cos^S 

tane  -  - * - 8  (6  -11) 

cosz8tan(8-6  ) 

8 

where  8  is  the  shock  angle  and  dg  the  shock  generator  angle. 

For  a  given  upstream  Mach  number,  this  equation  enables  the  variation  of 
e  with  6(  to  be  calculated. 
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McCabe  also  predicted  the  shock  generator  angle  to  induce  incipient  flow 
separation  for  given  upstream  Mach  number,  by  introducing  a  geometrical 
criterion  that  the  separation  appears  when  the  surface-flow  becomes  aligned 
with  the  shock,  i.e.  e  -  3  -  6g.  Using  (6  -11),  the  shock  generator  angle 
for  the  incipient  separation  is  calculated  by 

cos2(3  -  6  )  -  cosS  (6  -12) 

s 

which  is  the  form  given  by  Korkegi  (Ref. 9). 

McCabe's  theory  shows  fairly  good  agreement  with  the  experimental  maximum 
surface-flow  deflection  angles,  when  the  shock  is  not  very  strong,  (Fig. 39 
and  Fig. 57  mentioned  later).  However,  the  primary  assumption  of  infinitesimal 
boundary- layer  thickness  makes  it  impossible  to  include  any  details  of  the 
skewed  boundary  layer,  e.g.  the  velocity  profile.  Moreover,  his  criterion 
for  the  incipient  separation  cannot  always  agree  with  experimental  results, 
as  noted. 

For  the  skewed  boundary  layer,  a  more  realistic  (and  more  complicated) 
treatment  was  given  by  Lowrie  (Ref. 39  or  the  Review  by  Green,  Ref. 3). 

Lowrie's  analysis  of  the  cross  flow  in  the  boundary  layer  is  an  extension  to 
compressible  flows  of  the  suggestion  made  by  Johnston  (Ref. 40)  that  the  cross- 
flow  velocity  is  proportional  to  the  streamwise  velocity  defect  in  low- 
speed  boundary  layer.  Lowrie  employed  the  momentum  and  energy  equation  in 
his  analysis,  while  assuming  isentropic  flow  in  the  interaction  region,  and 
obtained  the  streamwise  and  cross-flow  velocities  in  an  outer  part  of  the 
boundary  layer  as  functions  of  the  pressure  rise  and  isobar  direction. 

Green  further  developed  this  theory  for  the  cross-flow  in  the  inner  part  of 
the  boundary  layer,  through  Lowrie's  'junction  point  criterion'  for  the 
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junction  between  the  outer  and  inner  layers.  Although  this  analysis 
gives  a  useful  theoretical  basis  for  the  cross  flow  behaviour*  it  includes 
a  number  of  physical  factors  which  are  very  difficult  to  either  estimate 
or  measure.  This  makes  comparison  with  experimental  data  very  dubious. 


6.2.3.  Angular-Momentum  Conservation  Concept 

The  analyses  of  McCabe  and  Lowrie  suggest  a  simple  approach  to  the  determination 
of  the  maximum  permissible  surface-flow  deflection  angle.  This  approach 
is  to  assume  that  angular  momentum  is  conserved  along  each  stream  tube  in  the 
side-wall  boundary  layer.  The  angular-momentum  conservation  concept  is  an 
extension  of  McCabe's  theory,  i.e.  replacement  of  the  total-circulation  in 
the  boundary  layer  by  the  angular  momentum  along  a  stream  tube.  This  change 
means  that  the  boundary  layer  condition  of  no-slip  at  the  wall  can  be 
satisfied.  The  viscous  flow  can  now  satisfy  Johnston's  cross-flow  model 
near  the  wall. 

Fig. 55a  shows  a  rectangular  stream  tube  altering  its  cross  section  and 
resultant  velocity  vector  whilst  approaching  a  surface-flow  convergence  line. 
(Although  the  stream  tube  is  drawn  to  run  along  the  wall  for  convenience, 
similar  stream  tubes  are  obtained  everywhere  through  the  skewed  boundary 
layer.)  The  stream  tube  is  defined  by  constant  mass-flux,  m  ,  through  it,  i.e. 

.  * 

m  ■  pU  hn  -  constant  (6  -13) 

where  p  is  density  and  U*  the  resultant  velocity,  while  h  and  n  are 
height  and  width  of  the  stream  tube  respectively.  The  alteration  of  the 
resultant  velocity  vector  brings  about  a  change  of  the  velocity  components 
parallel  and  normal  to  the  free  stream,  U,  V  and  W  respectively 
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corresponding  to  co-ordinate  axes  X,  Y  and  Z  (Fig. 55b).  If  it  is 
assumed  that  the  principal  angular-momentum  (H*)  related  to  the  resultant 
velocity  vector  (U*)  is  conserved  along  the  stream  tube  and  W  is  small 
enough  to  be  neglected  in  comparison  with  U  and  V  ,  then  U  and  V  are 
uniquely  related,  see  later  (equation  6  -15a).  The  principal  angular- 
momentum  flux  is  defined  using  the  mass  flux  and  the  vorticity,  (2  ,  as 


*  2. 

H  *  k  m  ft  ■  constant  , 

where  k  is  the  radius  of  gyration,  assumed  constant. 
Hence 


H 


.2. 

k  m 


3U 

3Z 


,  2. ,3U 
km(3Z 


Sing 


3V  r \ 
-  ^  cosg) 


tang 


3U  /3V 
3Z  /  3Z 


(6  -14) 


constant,  (6  -15a) 
(6  -15b) 


where  Z  is  the  displacement  from  the  wall  and  g  is  the  angle  of 
the  principal  angular  momentum  axis  to  the  free  stream.  Since  the  change 
of  the  streamwise  velocity  (U)  must  have  a  limit  depending  on  the  flow 
field  conditions  (such  as  pressure  rise),  so  the  cross-flow  velocity  (V) 
also  has  a  maximum.  From  these  limits  the  maximum  permissible  surface- 
flow  deflection  can  be  determined. 


For  the  analysis  presented  here,  it  is  more  convenient  to  use  a  co-ordinate 
system  which  is  valid  through  the  skewed  boundary  layer  over  the  whole 
interaction  region  and  which  needs  no  transformation  of  X  and  Y 
co-ordinates.  A  unit  flow  segment  of  the  stream  tube  is  defined  by 
satisfying  the  condition  that  its  area  projected  onto  the  wall  is  constant 
through  the  boundary  layer.  Equating  the  projected  area  to  that  at  the 
external  flow,  means  that 
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nU*At  -  n  U  At  (6  -16) 

e  e  . 

where  n  is  the  width  of  the  streamtube,  At  is  unit  time  and  suffix  e 

denotes  the  value  at  the  external  flow.  Substituting  (6  -16)  into 

(6  -13),  gives 

—  •  t-~  m  constant  (6  -17) 

P  n 
e  e 

along  the  stream  tube.  The  inviscid  external  flow  is  regarded  as 
parallel  to  the  wall,  so  that  an  inviscid  stream  tube  is  constant  in 
thickness  over  the  interaction  region.  Taking  hg  to  be  constant  in  this 
region  a  transformation  of  Z  is  introduced  from  (6  -17)  as  follows, 

z  -  (  (P/P  JdZ  (6  -18) 

'o 


6.2.4.  Analysis  of  the  Angular-Momentum  Conservation  Concept 
The  following  assumptions  are  made  in  this  analysis. 

(1)  The  boundary  layer  is  thin  and  the  flow  deflection  is  small,  so 
that  the  velocity  and  angular-momentum  components  normal  to  the 
wall  may  be  neglected  in  comparison  with  the  components  parallel 
to  the  wall. 

(2)  The  maximum  permissible  surface-flow  deflection  angle  is  determined 
by  the  greatest  change  of  flow  conditions  along  the  stress  tube, 
which  occurs  at  or  near  the  shock  location.  Therefore,  a  stream 
tube  crossing  the  shock  is  employed  in  the  analysis. 

(The  discussion  ought  strictly  to  be  made  upstream  of  the  surface- 
flow  convergence  line,  Fig. 55b). 

The  co-ordinate  axes  for  this  analysis,  x,  y  and  z,  together  with  the 
corresponding  velocity  components,  u,  v  and  w,  are  defined  in  directions 
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parallel  and  normal  to  the  external  flow  (Fig. 56b).  The  angular  momentum 
divided  by  mk^,  (H)  is  given  by  (6  -15),  as 

P-.  3u, 

Hi  ■  r  it  •  <6  -19a> 

le 

P7  3u,  3v, 

H2  -  pf  <-3T  8inC  -  "3z  cosC>  *  <6  -19b> 

2e 


tant; 


dz  /  9 z 


(6  -19c) 


where  suffices  1,2  and  e  denote  conditions  upstream  and  downstream 
of  the  shock  and  at  the  external  flow  respectively,  and  €  is  an  angle 
of  principal  angular-momentum  relative  to  the  external  flow,  while  z 
is  the  transformed  displacement  given  by  (6  -18) ,  as 

Z  "(  1^Pl/,ple^dzl  “(  2^P2^P2e^dz2  *  (6-20) 


Consider  an  idealized  boundary  layer  in  which  the  streamwise  velocity 
and  density  profiles  at  the  two  stations  are  similar  to  each  other, 

i,e’  u2  ■  Cuul  ’  (6  ~21a) 

P2  -  CpP1  (6  -21b) 

where  C  and  C  are  constants  determined  by  the  flow  field  conditions, 
up 

To  give  the  greatest  change  of  the  streamwise  velocity  profile  (see 
assumption  2) ,  these  constants  are  evaluated  from  the  oblique  shock 
relations 

Cu  -  cos8/cos(8  -  6s)  ,  (6  -22a) 

C  -  tan0/tan(0  -  6  )  (6  -22b) 

r  8 

where  6^  and  8  are  the  shock  generator  angle  and  shock  wave  angles 
respectively.  For  the  conservation  of  angular  momentum  Hj  -  H2  . 
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The  cross-flow  velocity  profile  is  obtained  in  the  following  manner. 
Substituting  (6  -21)  into  (6  -19)  gives 


sin£ 

3V2 

3z 


i  1  u  3z 


(6  -23a) 

(6  -23b) 


where  £  is  independent  of  z  through  the  boundary  layer. 

Since  v  «*  0  in  the  external  flow,  the  integration  of  equation  (6  -23b)  can 
be  written  as 


F7^ 


("l.‘  ul> 


(6  -24) 


which  is  the  cross-flow  velocity  profile.  However,  it  does  not 
give  zero  velocity  at  the  wall.  So,  it  is  a  solution  in  an  'outer 
layer'  only. 


According  to  Johnston's  cross-flow  model  in  the  'inner  layer'  near  the 
wall,  the  8treamwi8e  and  cross-flow  velocity  profiles  may  be  related  to 
each  other  by 

v 2  "  Auj  (6  ~25) 

where  A  is  independent  of  z  .  Substituting  (6  -25)  and  (6  -21)  into 
(6  -19)  gives 

tan?  ■  -  1/A  ,  (6  -26a) 

and  hence 


1  +  A 


(6  -26b) 


where  £  must  be  chosen  between  tt/2  and  ir  to  satisfy  equation  (6  -26). 
The  integration  of  equation  (6  -26)  gives 

u2  -  u1/  / 1  +  A2  . 


(6  -27) 
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To  satisfy  equations  (6  -27)  and  (6  -21a) , 
by 


A 


C  and  A  must  be  related 
u 


(6  -28) 


Hence,  the  inner  cross-flow  velocity  profile  is  given  by 


v 


2 


(6  -29) 


The  equations  for  the  inner  and  outer  profiles,  (6  -29)  and  (6  -24),  can 
be  matched  at  the  point  where 

u^  -  ulg/2  (6  -30) 

or,  using  (6  -21a) 

u2  *  u2e^2  *  <6  ~31^ 


The  direction  of  surface-flow  must  coincide  with  that  of  the  corresponding 
skin-friction,  so  that 


tane 


lim 

z-K> 


(III  1^1 

\  3z  /  3z 


(6  -32) 


where  the  surface-flow  angle,  e  ,  is  defined  relative  to  the  external 
flow.  Substituting  the  inn-layer  velocity  profiles  (6  -27)  and  (6  -29) 
into  (6  -32) ,  the  surface-flow  deflection  angle  is  given  by 


tane 


(8  - 
27 

cos  8 


V 


-  1  . 


(6  -33) 


For  a  given  upstream  Mach  number  (6  -33)  enables  the  variation  of  e 
with  6  to  be  calculated.  The  calculated  results  of  (6  -33)  for  various 
Mach  w bars  and  shock  generator  angles  are  presented  in  Fig. 57.  They 
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are  compared  with  experimental  surface-flow  convergence-line  angles 
obtained  from  both  the  present  measurements  and  some  earlier  work.  The 
experimental  angles  for  both  complete  and  incomplete  convergences  are 
indicated  using  solid  and  open  symbols  respectively.  The  figure  shows  that 
the  angle  predicted  by  (6  -33)  agrees  reasonably  well  with  the  condition 
for  appearance  of  the  complete  surface-flow  convergence.  In  the  figure, 
the  complete  convergence  line  angle  does  not  coincide  with  the  shock  wave 
angle . 


This  analysis  of  the  angular  momentum  conservation  concept  brings  out  some 
other  interesting  points.  Firstly,  the  location  of  the  junction  between 
the  inner  and  outer  parts  of  the  boundary  layer  occurs  at  a  point  where 
the  streamwi8e  velocity  (u2>  is  half  the  external-flow  velocity  (u2c) , 
as  given  by  (6  -31).  Provided  that  the  streamwise  velocity  profile  is 
conventional  and  decreases  monotonically  towards  the  wall  surface  then 
(6  -24)  and  (6  -29)  give  the  maximum  cross-flow  velocity  at  the  junction 
point  as 


(v2/u 


2e^max 


1  |cos2(B-62) 

2  /  Ta 

J  cos  B 


(6  -34) 


Experimental  data  taken  from  the  boundary-layer  surveys  of  Lowrie  (Ref. 39) 
and  Peake  (Ref. 13)  suggest  that  this  analytical  result  is  physically 
reasonable  as  shown  in  the  following  table: 
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REFERENCE 

TEST 

EXPERIMENTAL 

CORRESPONDING 

MAXIMUM 

CONDITIONS 

MAXIMUM 

STREAMWISE 

CROSS-FLOW  VELOCITY 

CROSS-FLOW 

VELOCITY  RATIO 

RATIO  FROM  EQUATION 

VELOCITY  RATIO 

(6  -34) 

Lowrie 

M  - 

2.5, 

0.15  -  0.18 

0.55  -  0.62 

0.16 

(Ref. 39) 

&  - 

s 

6° 

Peake 

M  - 

2, 

0.14  -  0.20 

0.45  -  0.52 

0.22 

(Ref. 14) 

6  - 

s 

8° 

M  - 

4, 

0.15  -  0.17 

0.45  -  0.50 

0.15 

<5  - 

s 

8° 

When  the  surface-flow  becomes  aligned  with  the  shock,  e  ■  8  -  <Sg,  so  that 


from  equation  (6  -33) 


2 


cos  e 


COSfJ 


which  is  the  same  result  as  obtained  by  McCabe. 

This  means  that  if  McCabe's  criterion  for  separation  is  accepted  then 
both  theories  predict  the  same  wedge  angle  to  cause  separation.  However, 
as  pointed  out  above,  McCabe's  criterion  is  not  accepted  here  and  a  better 
alternative  based  on  the  conservation  of  H*  is  shown  in  Fig. 57. 
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7. 


FLOW  FIELD  STRUCTURE  IN  THE  INTERACTION  REGION 

7.1  Flow  Structure  near  the  Shock  Generator 


The  oil-flow  pictures  taken  on  the  shock  generator  surface  show  that  two 
groups  of  compression  waves  affect  that  surface  and  expansion  waves  between 
them  accelerate  the  surface  flow  toward  the  shock  generator  root  (Figs. 37 
and  40) .  A  flow  mechanism  which  could  generate  these  compression  and 
expansion  waves  is  illustrated  in  Fig. 58a  which  shows  a  shock  envelope  around 
the  shock  generator,  together  with  a  typical  vapour-screen  picture  indicating 
the  sectional  view  of  the  flow  field  normal  to  the  free  stream.  The  side- 
wall  boundary  layer  forces  the  oblique  shock  to  become  detached  near  the 
shock-generator  leading-edge  root.  The  vapour-screen  picture  suggests  that 
the  shock  envelope  approaches  the  side  wall  in  the  form  of  a  lambda  shock. 


A  sectional  view  of  this  shock  envelope  in  the  x  -  z  plane  is  imagined  to 
be  as  drawn  in  Fig. 58b,  where  a  typical  oil-flow  picture  on  the  shock 
generator  surface  is  also  presented.  In  the  figure,  the  front  shock  originates 
from  a  group  of  compression  waves  springing  from  the  growth  of  the 
displacement  thickness  of  the  side-wall  boundary  layer.  This  front  shock 
results  from  adverse  pressure  propagation  through  the  subsonic  portion  of 
the  boundary  layer.  It  is  well  known  in  two-dimensional  flow  that  a  similar 
process  to  this  induces  a  very  thick  subsonic  layer  on  the  wall  surface. 
Although  the  present  flow  is  three-dimensional,  the  flow  in  the  x  -  z  plane 
which  cuts  the  shock  envelope  roughly  at  right  angles  seems  to  behave  in  a 
similar  way.  It  follows  that  a  recovery  in  thickness  of  the  side-wall 
boundary  layer  thins  the  thickened  subsonic  layer  while  shedding  expansion 
waves.  This  recovering  side-wall  boundary  layer  thins  locally  near  the  shock- 
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generator  leading  edge  root  (Fig. 58b).  The  second  group  of  compression  waves 
are  produced  by  the  corner  boundary  layer  growing  along  the  shock  generator 
root. 

As  indicated  in  the  oil-flow  picture  of  Fig. 58b,  the  first  and  second  groups 
of  compression  waves  form  the  outer  and  inner  compression  zones  on  the  shock 
generator  surface,  while  the  expansion  waves  give  rise  to  the  accelerated 
region  in  between. 

The  accelerated  region  drives  the  shock  generator  surface  flow  to  deflect 
toward  the  shock  generator  root.  Meanwhile,  the  growth  of  the  corner  boundary 
layer  with  a  thick  subsonic  region  pushes  some  shock-generator  surface  flow 
away  from  the  corner.  As  a  result  of  these  actions,  a  surface-flow 
convergence  (convergence  B,  Fig. 37)  is  formed  on  the  shock  generator  surface 
near  the  root. 


7.2.  'Double  Viscous-Layer  Model' 

Many  flow  field  models  of  the  glancing  shock  boundary-layer  interaction 
region  have  been  proposed  and  discussed  over  the  past  fifteen  years.  The  most 
popular  is  the  'vortex  dominated  model'  which  may  be  developed  from  an  idea 
based  on  a  conical  flow  field  model,  (see  Fig. 10).  This  model  consists  of 
a  'separation  vortex  (or  bubble)'  situated  between  a  'separation  line'  running 
a  little  upstream  of  the  shock  and  a  'reattachment  line*  lying  near  the  shock 
generator  root.  When  the  shock  generator  angle  is  fairly  large,  this  model 
could  be  suitable  for  interpreting  some  well-known  experimental  results. 


especially  the  surface  flow  patterns  which  were  most  often  available.  The 
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typical  interpretation  may  be  summarised  as  follows. 

(i)  A  side-wall  surface-flow  convergence  (convergence  A)  observed  in  the  oil- 
flow  pattern  a  little  upstream  of  the  shock  indicates  a  three-dimensional 
separation  line.  In  the  vicinity  of  this  flow  convergence,  the 
corresponding  surface  pressure  distribution  tends  to  have  a  small  peak  or 
flattened  part  which  suggests  a  rapidly  thickned  separated  layer  being 
unable  to  withstand  a  strong  adverse  pressure  gradient. 

(ii)  Oil  streak  lines  diverging  near  the  shock  generator  root  are  caused  by  a 
reattachment  line  just  behind  the  'separation  vortex'.  Near  the  re¬ 
attachment  line,  the  impingement  of  air  entrained  by  the  vortex  yields 

a  peak  in  the  heat  transfer  distribution  and  a  crest  in  the  surface 
static  pressure  plateau. 

However,  this  flow  field  model  does  not  explain  all  the  present  experimental 
results  obtained  in  the  two  separate  test  programmes,  as  mentioned  below. 

Oil-Flow  Pictures 

The  'complete  surface-flow  convergence  (convergence  A)'  on  the  side  wall  may  be 
taken  as  evidence  of  'separation'  appearing  in  this  region.  However,  before 
the  appearance  of  the  separation  line  on  the  side  wall,  there  is  already  evidence 
of  a  small  separation  vortex  near  the  shock  generator  root,  as  indicated  by  the 
typical  oil-flow  pattern  of  flow  reattachment,  (Fig. 37  and  Section  4.3.1.). 


Vapour-Screen  Photographs 

The  vapour-screen  photographs  show  that  a  small  vortex  is  swirling  near  the 
■hock  generator  root  and  a  thick  shadowy  zone  is  located  between  the  shock 


-  81  - 


generator  root  and  the  side-wall  surface-flow  convergence  line  (convergence  A) . 
However,  these  features  are  observed  at  all  shock  generator  angles  so  that 
no  particular  change  of  pattern  suggesting  the  appearance  of  a  'separation 
vortex'  can  be  detected.  (Figs. Ale  -  Ale  and  Section  A. 3. 3.). 

Heat  Transfer 

Even  before  the  'complete  surface-flow  convergence'  (convergence  A)  takes 
place  on  the  side  wall,  the  heat  transfer  distribution  has  two  'peaks' 

(at  least),  ahead  of  the  shock  and  near  the  shock  generator  root.  (Figs. 

AAe  -  AAg  and  Section  A. A. 2.). 

A  flow  field  model  which  is  consistent  with  these  flow  features  is  the 
'double  viscous-layer  model'  illustrated  in  Fig. 59  and  discussed  below. 

(1)  The  Viscous  Flow  Separating  from  the  Shock  Generator  Surface  and 
Reattaching  to  the  Side  Wall  Surface 

The  oil-flow  pictures  and  the  vapour-screen  photographs  for  the  shock 

generator  surface  at  any  angle  to  the  free  stream,  show  that  a  spanwise 

surface  flow  on  the  shock  generator  is  accelerated  toward  the  side  wall 

through  the  'accelerated  region'  (Fig. 37  or  AO)  and  a  small  vortex  is 

swirling  near  the  shock  generator  root  (Fig.Al).  It  is  considered  that  the 

small  vortex  is  caused  from  the  spanwise  surface  flow  separating  at  its 

convergence  line  (convergence  line  B)  from  the  shock  generator  surface  and 

reattaching  to  the  side  wall  at  the  divergence  line  near  the  shock  generator 

root  (Fig. 59).  This  flow  structure  indicates  that  some  viscous  flow  growing 

on  the  shock  generator  surface  comes  into  the  side  wall  surface,  by  way  of 

the  edge  of  the  small  'roller  vortex'  in  the  corner.  This  impingement  of 

the  viscous  flow  produces  a  'close  packed'  surface-flow  on  the  side  wall 
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near  reattachment,  as  suggested  by  the  oil  streak  lines  with  fine  intervals 
in  this  region.  The  impingement  gives  rise  to  the  heat  transfer  peak 
(peak  (D  near  the  shock  generator  root  (Figs.44e  -  44g) .  This  peak  is 
therefore  not  connected  with  the  main  separation  (at  the  convergence  line  A) , 
namely  the  side-wall  boundary-layer  separation. 

(2)  The  Viscous  Layer  Induced  on  the  Side  Uall  Surface 

On  the  side-wall  surface,  the  impingement  of  the  viscous  flow  from  the  shock 
generator  surface  produces  a  viscous  layer  in  the  vicinity  of  the  reattachment. 
This  'induced  layer'  is  accelerated  toward  the  shock  by  the  surface  static 
pressure  gradient  (Fig. 43)  and  crosses  the  path  of  the  side-wall  boundary 
layer  coming  into  the  region,  so  drawing  a  boundary  line  on  the  side  wall 
surface  between  itself  and  the  side-wall  boundary  layer  as  illustrated  in 
Fig. 59.  This  process  causes  the  surface-flow  convergence  (convergence  A) 
near  the  shock  location  (Fig. 37). 

(3)  The  Side-Wall  Boundary  Layer  upstream  of  the  Inviscid  Shock 
Upstream  of  the  calculated  shock,  the  surface  static  pressure  distribution 
shows  an  extensive  region  of  adverse  pressure  gradient  propagated  through 
the  subsonic  region  of  the  side  wall  boundary  layer  (Fig. 43).  This  pressure 
propagation  causes  the  shock  to  take  a  lambda  configuration  near  the  side 
wall  surface,  as  shown  in  the  vapour-screen  pictures  (Fig. 41).  As  a  result 
the  gradual  compression  of  the  side-wall  boundary-layer  flow  tends  to  increase 
the  heat  transfer.  As  the  convergence  line  A  is  approached  so  the  surface 
streamlines  swing  round  and  pack  closer  together  which  gives  rise  to  the 

heat  transfer  peak  (A)  . 


(4)  Boundary  Region  between  the  Side-Wall  Boundary  Layer  and  the  Induced  Layer 

When  the  shock  generator  angle  is  not  large,  the  surface  static  pressure 
gradient  on  the  side  wall  is  not  very  steep  (Figs. 43a  -  43b  and  43d  -  43e). 

The  induced  layer  cannot  develop  fully  over  the  interaction  region.  The 
side-wall  boundary  layer  may  turn  around  the  edge  of  the  induced  layer  while 
keeping  its  surface  stream  lines  in  contact  with  the  side  wall  surface 
(Figs. 59  and  60a).  This  process  yields  an  attached  flow  field  showing  the 
incomplete  oil-flow  convergence  (convergence  A)  on  the  side  wall  surface  near 
the  shock  location  (Figs. 37a  and  37d). 

With  increasing  shock  generator  angle,  the  steepened  surface  static  pressure 
gradient  (Figs. 43c  and  43f)  accelerates  the  induced  layer  further  upstream 
so  that  the  front  edge  of  the  induced  layer  presents  a  large  angle  to  the  free 
stream.  In  this  case,  the  side  wall  boundary  layer  can  ride  over  the  induced 
layer  and  its  surface  stream  lines  separate  from  the  side  wall  surface 
(Fig8.59  and  60b).  (Refer  to  Section  6.1).  This  situation  represents  a  flow 
separation  similar  to  Maskell's  'ordinary  separation',  and  the  corresponding 
oil-flow  pictures  show  the  complete  surface-flow  convergence  (convergence  A) 
on  the  side  wall  upstream  of  the  shock  (Figs. 37b  -  37c  and  37e  -  37f ) . 

When  this  ordinary  separation  takes  place,  the  side-wall  boundary  layer  and 
some  viscous  flow  of  the  induced  layer  are  transformed  into  a  separated  free 
layer  (Fig. 60b).  This  process  looks  similar  to  'separation  in  the  vortex 
dominated  model'.  However,  the  present  experimental  results  show  a  different 
separation  pattern  from  that  of  the  vortex  dominated  model  in  two  ways. 

First,  the  appearance  of  the  surface-flow  divergence  pattern  (showing  a 
flow-reattachment)  near  the  shock  generator  root,  as  well  as  the  heat  transfer 
peak  (B)  near  it,  has  no  connection  with  occurence  of  the  main  separation 
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(Che  complete  convergence  A)  on  the  side  wall,  as  noted  in  Figs. 37  and  44. 
Secondly,  the  present  vapour- screen  pictures  do  not  show  any  particular  flow 
features  to  give  evidence  of  the  separation  vortex  originating  from  the 
main  separation  on  the  side  wall.  Fig. 41.  Thus,  it  is  considered  that  the 
side-wall  boundary  layer  separation  does  not  produce  a  strong  vortex  but  a 
vortical  free  layer,  and  the  free  layer  near  the  shock  generator  can  become 
involved  with  the  small  vortex  swirling  at  the  shock  generator  root  (Fig. 60b). 


7.3.  An  Hypothesis  for  the  Appearance  of  Flow  Separation 

In  the  'double  viscous-layer  flow  field  model'  proposed  in  the  last  Section, 
the  occurrence  of  side-wall  boundary-layer  separation  is  not  directly 
connected  with  the  surface-flow  angle  on  the  side  wall  exceeding  the  shock 
angle.  The  discussion  here  has  been  based  on  experimental  results  showing 
'conical'  flow  features;  e.g.  the  interaction  region  grows  in  a  radial 
direction.  The  fact  that  the  present  interaction  region  is  'conical'  in 
nature  casts  doubt  on  McCabe's  criterion  for  incipient  separation  which  was 
based  on  a  quasi- two-dimensional  view  of  the  flow. 

This  point  may  be  confirmed  by  Fig. 61  in  which  a  correlation  of  the  surface 
convergence- line  angle  with  the  appearance  of  side-wall  boundary-layer 
separation  is  presented,  for  both  the  experimental  data  of  the  present  study 
and  a  number  of  other  references.  In  the  figure,  the  ratio  of  experimental 
convergence- line  angle  to  the  theoretical  prediction  (the  angular  momentum 
conservation  concept  proposed  here  and  McCabe's  theory)  is  plotted  versus 
the  shock  strength  represented  by  static  pressure  ratio  across  the  shock. 


Following  McCabe's  criterion,  Korkegi  proposed  a  pressure  ratio  of  1.5 
to  induce  incipient  separation.  However,  in  the  current  work  the  incomplete 
type  of  surface-flow  convergence  can  appear  even  when  the  pressure  ratio 
is  far  beyond  1.5.  The  condition  for  the  appearance  of  separation, 
defined  by  complete  surface-flow  convergence  on  the  side  wall,  occurs 
when  the  surface-flow  angle  on  the  convergence  line  exceeds  a  maximum 
permissible  angle.  This  angle  is  better  predicted  by  the  angular  momentum 
conservation  concept,  than  by  McCabe's  theory,  though  neither  theory  gives 
particularly  good  agreement  with  the  experimental  data,  see  Fig. 61. 

The  following  hypothesis  is  proposed.  An  ordinary  separation  occurs  on 
the  side  wall  in  the  interaction  region,  when  the  induced  layer  forces 
the  surface  stream  lines  of  the  side-wall  boundary  layer  to  deflect 
beyond  a  maximum  permissible  angle.  This  angle  is  not  the  shock  angle. 
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CONCLUSIONS  AND  REMARKS 

(1)  The  three  dimensional  glancing  shock  turbulent  boundary-layer  interaction 
region  consists  of  two  different  viscous  layers;  i.e.  (i)  the  side-wall 
boundary  layer  growing  along  the  flat  wall  and  (ii)  the  induced  layer 
originating  on  the  shock  generator  surface  near  the  root  and  crossing  the  path 
of  the  side-wall  boundary  layer. 

(2)  As  long  as  the  surface  stream  lines  of  the  side-wall  boundary  layer 
are  pliable  enough  to  be  bent  along  the  edge  of  the  induced  layer,  no  flow 
separation  appears  on  the  side  wall  even  if  the  surface-flow  deflection 
exceeds  the  shock  angle.  An  ordinary  separation  takes  place  when  the  induced 
layer  forces  the  side-wall  surface  stream  lines  to  deflect  beyond  a  maximum 
permissible  angle. 

(3)  McCabe' 8  criterion  is  not  suitable  for  deciding  the  incipient  flow 
separation  on  the  side  wall,  since  the  interaction  is  not  quasi-two-dimensional. 
As  defined  here  separation  occurs  when  the  surface-flow  deflection  exceeds 

the  maximum  permissible  angle.  This  angle  may  be  determined  by  an  angular 
momentum  relation  for  the  side-wall  boundary  layer. 

(4)  If  the  heat  transfer  distribution  in  the  interaction  region  is  to  be 

measured  under  conditions  when  the  temperature  difference  (Tw  -  Tftw)  is  small, 

it  is  essential  to  take  into  account  the  local  variations  in  both  T  and  T  , 

w  aw 

particularly  in  the  neighbourhood  of  the  shock. 
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Figure  1.  An  example  of  high  local  heating. 
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Figure  2.  An  example  of  shock-induced  boundary 

layer  interactions. 
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Figure  3.  Typical  two  dimensional  shock 
boundary -layer  interactions. 
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Figure  4.  Some  experimental  and  theoretical 
results  of  two  dimensional  shock  boundary- 

layer  interaction. 
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Figure  6.  Masked ‘s  separation  models. 
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Figure  7.  Quasi -two-dimensional  flow  field. 
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Figure  7.  Quasi -two -dimensional  flow  field 
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Figure  15.  Model  geometry  in  the  2-5  x  2  5  inch  tunnel 
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Figure  17.  Roughness  elements  to  modify  the 

test  boundary  layer. 


Figure  18.  Working  section  of  the  9x9 inch  tunnel. 


Figure  23.  Arrangement  of  surface  static  holes. 
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Figure  24.  Temperature -time  histories  for  the 
slug  calorimeter  technique. 
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Figure  29.  Temperature -time  history  of  a  thin -film 
gauge  in  the  continuous  running  tunnel. 
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Figure  30.  Thin -film -gauge  calorimeter  installation. 
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(b)  Electronics  installation. 


Figure  30.  Thin-film-gauge  electronics  installation. 
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Figure  33.  Encapsulated  liquid  crystal  technique. 
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Figure  35.  Co-ordinate  systems  for  experimental 

data  presentation. 
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Fig.  37  Contd. 
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Fig.40  Detailed  Oil-Flow  Patterns  on  the  Shock 
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Fig.  40  Contd. 
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Fig.  41  Contd. 
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Figure  42.  Shock  locations  as  measured  from 
the  schlieren  and  vapour-screen  pictures. 


Figure  43.  Isobars  on  the  side -wall  surface. 


Figure  43.  Isobars  on  the  side-wall  surface 
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Figure  44.  Heat  transfer  on  the  side- wall  surface. 


Figure  44.  Heat  transfer  on  the  side-wall  surface 
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Figure  48.  Static-pressure  and  heat -transfer 

distribution. 
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Figure  49.  Correlation  between  static  pressure 

and  heat  transfer. 
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Figure  49.  Correlation  between  static  pressure 

and  heat  transfer. 
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Figure  49.  Correlation  between  static  pressure 

and  heat  transfer. 
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Figure  49.  Correlation  between  static  pressure 

and  heat  transfer. 
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Figure  50.  Viscous -layer  survey  results. 
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Figure  53.  Typical  heat  transfer  distributions. 


i 


•r  4*  '«•> 


(a)  McCabe's  Theory 
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Figure  56.  Co-ordinate  systems  for  surface 

flow  deflection  analyses. 
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Figure  56.  Row  structure  near  the  shock 

generator  root. 
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Experimental  Angle  Experimental  Angle 

McCabe's  Prediction  The  Present  Theory 


1.5 

1.0 

0.5 

1.5 

1.0 

1.0 

12  3  4 

Shock  Strength  (Pressure  Ratio  across  the  Shock) 

o  cf  Model  A  and  B 
a  d  Ref.  11  (Model  1  and  2) 

0  4  Ref.  14  (Ms  2  and  4) 


open  symbol:  incomplete  conv. 
solid  symbol:  complete  conv. 


Figure  61.  Correlation  of  convergence  line  angle 
with  the  appearance  of  flow  separation 


APPENDIX  A 


THE  DERIVATION  OF  THE  BOUNDARY  LAYER  PROFILES 


The  boundary  layer  was  surveyed  using  Pitot  and  static-pressure  probes. 

These  two  pressures  (and  the  measured  stagnation  temperation,  T  )  gave  local 
Mach  number,  density  and  velocity  profiles  in  the  boundary  layer  using 
the  relations  outlined  below. 


The  supersonic  boundary- layer  flow  cannot  be  regarded  as  isentropic  and  the 


measured  local  Pitot  pressure,  pp  ,  is  the  stagnation  pressure  downstream 
of  a  normal  shock  induced  by  the  probe  head.  Hence  the  local  Mach  number, 
M  ,  must  be  derived  from  Rayleigh's  supersonic  Pitot  formula,  i.e. 


JL 

P„ 


/  J]L  M2  - 

V  y+1  Y+l / 


V  V+l 


(^H2) 


T*l' 

y/(y-i) 


(A  -  1) 


where  p  is  the  measured  local  static  pressure  and  Y  is  the  ratio  of 
specific  heats  assumed  to  be  1.4  for  air. 

The  calculated  Mach  number  profile  can  then  be  used  to  derive  the  density 
profile  p(y)  from  the  equation  of  state, 

P  -  p/RT  ,  (A  -  2) 

where  R  is  the  gas  constant  (287  J/kg.°K)and  T  is  the  local 
temperature  calculated  from 


T/T  -  1/(1  +  ^  M2) 

O  i 


(A  -  3) 


Finally,  the  resultant  velocity  profile,  U  ,  is  derived  from 

IpU*2  -  ^  M2  (A  -  4) 


which  is  the  definition  of  dynamic  pressure 


APPENDIX  B 


HEAT  TRANSFER  MEASUREMENTS 


B.l.  Quasi-Transient  Calorimeter  Techniques 

There  are  many  different  techniques  for  measuring  the  heat  transfer  between 
a  solid  wall  and  a  moving  airstream.  They  fall  roughly  into  three  different 
categories,  according  to  the  temporal  behaviour  of  the  measured  temperatures 
(e.g.  wall,  adiabatic-wall,  etc.)  which  are  to  be  reduced  to  the  heat 
transfer.  The  techniques  are 

(1)  the  transient  technique  in  which  all  the  temperatures  are  variables  of 
time, 

(2)  the  quasi-transient  technique  in  which  only  the  local  wall  temperature 
is  regarded  as  variable  with  time,  and 

(3)  the  steady-state  technique  in  which  the  temperatures  have  no  change 
in  time  but  in  space. 

Most  of  these  techniques  are  further  divided  into 

(1)  surface- thermometer  techniques,  and 

(2)  calorimeter  techniques. 

The  difference  between  the  surface-thermometer  and  calorimeter  techniques 
is  small.  However,  for  the  calorimeter  technique,  a  small  'calorimeter  area' 
must  be  locally  and  separately  situated  on  the  wall  surface. 

To  measure  the  heat  transfer  distributions  in  the  present  shock  boundary-layer 
interaction  region,  the  quasi-transient  calorimeter  technique  was  used  for 
both  test  programmes.  However,  the  details  of  the  technique  used  were 
different  between  the  test  programmes.  This  was  because  of  basic  differences 
in  the  test  conditions  between  them,  i.e.  intermittent  and  continuous  test 
conditions.  In  the  2.5  x  2.5  inch  tunnel  test  prograane,  a  modified 
'Princeton  University  calorimeter  technique'  (Appendix  B.2)  was  employed. 


But  for  the  9x9  inch  tunnel  programme  a  'thin-film-gauge  calorimeter 
technique'  (Appendix  B.3)  was  devised  and  developed.  The  general  principle 
of  these  calorimeter  techniques  is  as  follows. 


The  thermal  field  where  the  heat  transfers  between  the  solid  wall  and  the 
air  occurs  is  three-dimensional.  However,  this  thermal  field  may  be  regarded 
as  one-dimensional,  if  the  component  of  the  heat  flux  parallel  to  the 
wall  is  negligible  in  comparison  with  the  normal  component  (Fig.Bl). 

The  thermal  balance  at  the  wall  surface  may  be  written  using  the  heat-transfer 
rate  (q)  or  the  heat-transfer  coefficient  (ha> ,  as  follows: 


q  * 

in  the  'wall  domain'. 


T 

k  i-s 

w  3z 
and 


(B  -  1) 
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-  k 


3T 
_ a 

a  3z 


(B  -  2) 


or 

«  ■  V(V.-o  -  T.„»  <,-3> 

in  the  'air  domain',  where  T,  k  and  z  are  temperature,  thermal  conductivity 
and  distance  from  the  wall  surface  respectively,  while  suffices  w,  aw  and 
a  denote  conditions  of  the  wall,  the  adiabatic  wall  and  the  air.  In  the 
wall  domain,  the  one-dimensional  heat  equation  may  be  written  as 


32T 


3z 


i  3t 

1  w 


a  3t 
w 


where  t  is  time  and 


(B  -  4) 

is  thermal  diffusivity  of  the  wall  material 


defined  by  in  terms  of  thermal  conductivity,  k^  , 
specific  heat,  cw>  (<*w  -  k/Vw0 


density,  ptf  and 


Equation  (B  -  4)  can  be  integrated  with  respect  to  y  as 


3T 

w 

dz 


T  « 

_L _ m 

a  dt 
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(B  -  5) 


where  T  is  a  mean  temperature  in  a  small  domain,  Z,  and  L  is  defined  by 
n 


dz 


(B  -  6) 


Substituting  (B  -  5)  into  (B  -  1)  gives 


dT 

-0  c  L  — - 
w  w  dt 


(B  -  7) 


It  must  be  noted  that  p^c^L  is  an  equivalent  thermal-capacity  of  the  wall 
surface. 


B.2.  Slug  Calorimeter  Technique 

It  is  assumed  that  a  small  'slug*  is  situated  on  the  wall  surface  but 
thermally  insulated  from  it  (Fig.B2).  When  this  slug  is  used  as  a 
calorimeter  with  known  thermal  capacity,  C  ,  and  surface  area,  A  ,  equation 
(B  -  7)  may  be  rewritten  as 


where  the  suffix  s  denotes  conditions  of  the  slug.  Thus,  the  heat 
transfer  rate  from  the  slug  to  the  air  is  directly  obtained  by  measuring 
the  time  history  of  T  . 

B8 


From  (B  -  3)  it  follows  that 


when  ha  is  assumed  independent  of  temperature  and  the  thermal  resistance 

in  the  slug  is  small  enough  for  T  to  represent  the  surface  temperature 

ms 

of  the  slug.  Hence,  it  is  clear  from  (B  -  9}  that  the  heat-transfer  rate 
from  the  wall  to  the  air  (q)  can  be  determined  by  measuring  another  heat- 
transfer  rate  from  the  slug  (qg) ,  when  ■  (T^y-O,  because  under 

these  conditions  q^/q  ■  1. 

B.3.  Thin-Film-Gauge  Calorimeter  Technique 

Let  us  consider  a  thin-f ilm-gauge  installation  (Fig.B3),  made  up  from  a  thin 
Pyrex  glass  plate  (the  outer  wall  layer) ,  and  a  thick  aluminium  slab  (the 
inner-wall  material).  The  thin-f ilm  gauges  are  sputtered  on  the  glass- 
plate  surface.  This  configuration  can  provide  a  calorimeter  for  each  thin- 
film  gauge,  if  use  is  made  of  the  differences  in  thermal  capacity  and 
conductivity  between  the  inner  and  outer-wall  materials.  The  thermal- 
conductivities  of  Pyr ex-glass  and  aluminium  are  0.014  and  2.4  J/cm.sec.°K 
respectively.  The  much  larger  thermal-conductivity  and  capacity  of  the 
aluminium  wall  ensures  that  the  temperature  of  the  inner  wall  is  approximately 
uniform  and  constant.  If  the  glass-plate  surface  is  heated  locally  by 
the  thiu-film  gauge,  a  temperature  distribution  is  formed  in  the  glass-plate 
layer,  as  shown  in  Fig.B4.  This  small  heated  part  of  the  glass-plate  is  the 
calorimeter.  The  advantages  and  disadvantages  of  such  a  calorimeter,  as 
used  in  the  quasi-transient  technique,  were  described  in  Section  3.5.2. 

The  thermal  field  for  this  calorimeter  has  two  different  heat-transfer  systems 
(Fig.B3b),  i.e.  heat  transfer  from  the  calorimeter  region  of  the  wall  surface 


Co  the  air  (q^)  and  heat  transfer  from  the  inner  wall  to  the  calorimeter  (q^) 
Therefore  equation  (B  -  7)  may  be  written  as 

^  dT 


-  -  C 


dt 


(B  -10) 


where  is  the  gauge  temperature  and  C  is  an  equivalent  thermal- 
capacity  per  unit  area  of  the  calorimeter  (which  is  discussed  in  Appendix  B4) 
Equation  (B  -  3)  applied  locally  gives 


-  h  (T  -  T  ) 
a  g  aw 


(B  -11) 


h  (T 
w  wr 


(B  -12) 


where  h  and  h  are  heat-transfer  coefficients  corresponding  to 

8  W 

qa  and  qw  respectively  and  is  the  inner-wall  temperature  which  is 
used  as  a  reference  value. 


From  the  thermal  balance  at  the  wall  surface,  it  is  clear  that,  only 
when  dTg/dt  •  0,  does 

q  "  qa  "  %  (B  '13) 

where 

Tg  -  <Vz-0  (B’14) 

Hence  this  type  of  calorimeter  cannot  directly  measure  the  heat- 
transfer-rate  from  the  calorimeter  to  the  air,  but  it  can  be  used  to 
measure  the  heat-transfer-coefficient  in  the  following  manner. 

Substituting  (B  -11)  and  (B  -12)  into  equation  (B  -10)  gives 
dT 

c‘  it  *  *  v«,  -  v  -  w  -  *«>  • 


(B  -15) 


If  and  are  assumed  independent  of  temperature,  then  this  equation 

has  only  one  solution  which  satisfies  the  initial  condition  that 

T  •  T  .  when  t  *  0,  i.e. 
g  gi 


T  -  (T  )  . 

K  w  Z»C 


V  (Vz-< 


/  h  +  h  \ 

(-•Vs- *) 


(B  -16) 


where  (B  -14)  is  used  to  determine  the  final  value  of  T^,  i.e.  (T^^q 

which  is  a  function  of  T  . 

aw 

It  must  be  noted  that  the  heat-transfer-coefficient  appears  as  a  term  in 


the  time  constant  of  equation  (B  -15)  and  the  adiabatic  wall  temperature 


where  s  and  Z  are  the  integration  domains  defined  by  the  surface 
area  of  the  calorimeter,  and  its  depth,  respectively. 


*00 


Strictly,  discussion  of  the  transient  temperature-distribution  involves 
taking  account  of  the  thermal  penetration  depth,  2 /at  where  a  is  the 
thermal  diffusivity.  In  the  present  case,  the  penetration  depth  for 
Pyrex-glass  is  0.07  inches  for  a  heat  pulse  of  1  second  duration,  the  glass- 
plate  thickness  is  0.058  inches  and  the  time  constant  of  equation  (B  -16) 
is  of  the  order  of  1  -  10  seconds.  Hence  it  is  considered  that  ignoring 

* 

the  penetration  depth  does  not  produce  noticeable  errors  in  estimating  C  . 

It  can  be  assumed  that  the  temperature  in  the  inner  wall  is  uniform  and 
constant  (T^) .  Then,  under  nearly  steady  conditions  (dTw/dt  =  0) ,  the  one¬ 
dimensional  heat  equation  (B  -  4)  gives  a  linear  solution  for  the 
temperature  distribution  through  the  glass  plate,  i.e. 

T»  •  -IV  <Vz.ol z/d  <»-19> 

where  d  is  the  thickness  of  the  glass  plate.  This  temperature  distribution 
suggests  some  types  of  thermal  models  to  determine  the  equivalent  thermal 
capacity. 

The  simplest  model  is  that  the  thin-film  gauge  provides  constant  temperature 
over  its  surface  and  (B  -19)  is  applied  throughout  the  whole  region, 
(Fig.B4a).  In  this  case,  the  equivalent  thermal  capacity  is  given  by 

C*  -  |pyCwd  (B  -20) 

However,  this  value  of  the  thermal  capacity  gives  poor  agreement  with 
the  experimental  results  for  bare-plate  heat-transfer  measurements. 

Moreover,  it  is  physically  unrealistic  to  have  discontinuities  in  the 
temperature  distribution.  The  best  agreement  with  the  experimental  results 


for  a  bare  flat  plate  is  obtained  using  the  'parabolic  thermal-model', 
as  described  in  Section  3.5.2.  This  model  has  a  parabolic  surface-temperature 
distribution  across  the  width  of  the  gauge,  as  proposed  by  Busing  (Ref. 28), 
and  a  parabolic  region  of  thermal  influence  of  maximum  depth  d,  (Fig.B4b). 

The  equivalent  thermal  capacity  of  this  model  is 

C*  -  ^  p  C  d  (B  -21) 

When  the  glass-plate  surface  is  heated  by  the  thin-film  gauge  and  the 
temperature  distribution  in  the  calorimeter  reaches  a  steady  state,  the 
maximum  temperature  in  the  glass-plate  layer  is  at  the  surface.  The 
temperature  distribution  through  the  layer  is  given  by  (B  -19) .  When  the 
'heating  process'  is  changed  to  the  'cooling  process'  to  indicate  the 
gauge-temperature  time-history  corresponding  to  (B  -16) ,  the  temperature 
distribution  through  the  layer  changes  and  the  maximum  temperature  will 
occur  at  some  point  within  the  layer.  The  thin-film  gauge  on  the  wall 
surface  operates  with  this  'locally'  transitional  surface-temperature 
variation  and  indicates  a  'transitional  period'  in  its  temperature  time 
history,  see  Fig. 29a.  After  this  temperature  re-distribution  in  the 
calorimeter,  it  is  difficult  to  estimate  the  equivalent  thermal  capacity. 
However,  if  the  heated  region  of  the  glass  plate  is  so  thin  (compared 
with  the  thermal  penetration  depth)  that  the  calorimeter  region  can  be 
treated  as  a  lump,  then  it  can  be  assumed  that  the  equivalent  thermal 
capacity  estimated  initially  by  (B  -18)  is  still  valid  even  after  the 
re-distribution. 
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